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By W. Z. Cuten’, L. INFELD, J. R. POUNDER, A. F. STEVENSON, 
AND J. L. SYNGE 


Abstract 


Part I deals with the problem of determining the field due to a source of 
radiation inside a semi-infinite rectangular wave guide closed at one end by a plug, 
the current distribution in the source being regarded as known. Both the 
walls of the guide and the plug are treated as being perfectly conducting. Three 
different methods of solving the problem are given. The radiation resistance 
is then deduced from energy considerations. In particular, an expression for 
the radiation resistance of a linear antenna perpendicular to the wider face of 
the plug, fed at the point of entry, is derived, it being assumed that the antenna 
current is sinusoidal and that only the fundamental H-wave is transmitted by 
the guide. 


In Part II, one of the methods of paper I is extended to the case of a guide of 
arbitrary cross section, and the general problem of the calculation of radiation 
resistance and reactance is discussed. 


In Part III, a number of formulae for the radiation resistance of antennae of 
various shapes, with various assumed current distributions, in rectangular and 
circular guides, are given. 

In Part IV, explicit calculations for the impedance of a linear antenna in a 
rectangular wave guide are given. Further numerical calculations relating to the 
same problem, from the point of view of matching and sensitivity, have been 
made by Messrs. Chien and Pounder, but are not reproduced here. 


General Introduction 


The four parts of this paper are based on four of the reports which were 
issued by the Special Committee on Applied Mathematics of the National 
Research Council, the authors of which were, at the time, members of the 
Department of Applied Mathematics of the University of Toronto. The 
four parts deal with rather closely related topics in the field of wave guide 
theory which are all of a certain mathematical interest. 

They have been reproduced with only slight alterations from their original 
form, and no attempt has been made to refer to more recent work, although, 
so far as the authors are aware, most of the results contained herein have not 
hitherto been published. References are given at the end of each part 
separately. The figures and equations are also numbered separately for 
each part. 

The authors wish to acknowledge the help received from the Field Station 
of the National Research Council. 

1 Manuscript received October 29, 1948. 
Contribution from the Department of Mathematics, University of Toronto, Toronto, Ont. 


Formerly of the University of Toronto. 
Institute of Advanced Studies, Dublin. 
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PART I. RADIATION FROM A SOURCE INSIDE A PERFECTLY 
CONDUCTING WAVE GUIDE OF RECTANGULAR SECTION 


By L. INFELD, A. F. STEVENSON, AND J. L. SYNGE 


1. Introduction 


We consider a wave guide of rectangular section, extending to infinity in 
one direction and closed in the other direction by a plug. The guide and the 
plug are assumed to be perfect conductors. 

A source of energy is placed inside the guide. First we take this source 
to be a dipole, arbitrary in position and direction. We find the field produced 
inside the guide by the dipole, particularly the field a long way down the 
guide. This gives the mean flux of energy at infinity in the guide. Since the 
walls and plug are perfectly conducting, no energy can pass into them or out 
of them, and so this flux of energy is the power output of the dipole. Dividing 
this output by the mean value of the square of the current in the dipole, we 
obtain the radiation resistance. 

Once the field due to a dipole has been found, the field due to a finite 
antenna follows by integration. The power output of the antenna is 
obtained from the field a long way down the guide, and the radiation resist- 
ance of the antenna is given by dividing the output by the mean value of the 
square of the current entering the antenna. We thus find the resistance, but 
not the reactance. The Heaviside (or rational) system of units is used 
throughout, the final results being converted into ohms. 

We have succeeded in solving our problem by three independent methods, 
all of them giving the same result. Since we do not know which of these 
methods will prove most useful in the calculation of reactance, we give all 
three of them here. The first two methods are included in the main body of 
the paper, the third in Appendix B. 

These three different methods can be characterized as follows. 

The first two methods are based on the idea of images. The first (Section 3) 
is the more intuitive one; it treats the dipoles as singularities throughout, 
and gives the field only at great distances from the lattice. The second 
(Section 4) is the more rigorous one; it gives the damped waves near the lattice 
as well as the field at great distances; it is characterized by a smoothing process 
which enables us to replace summations by integrations. Finally, the third 
method (Appendix B) avoids the use of images altogether (except for images 
in the end of the guide) and calculates the field not by means of the Hertzian 
vector, but directly from Maxwell’s equations, assuming a known current 
distribution. 

Our method is sufficiently general to cover any rectangular section for the 
guide and any form of antenna, assumed to be a very thin wire. Details are 
developed only for a guide which cuts off all radiation except the Ho wave; 
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for this case a general formula for radiation resistance is given, the antenna 
remaining arbitrary. In particular we take the antenna to be straight and 
perpendicular to the wider face of the guide, being fed at the point of entry. 
Then the radiation resistance is 





(1.1) = “0 a ‘os? ~~ sin? ory tan? A ohms, 
where a = _ longer edge of cross section in centimeters, 
b = shorter edge of cross section in centimeters, 
AX = wave length in centimeters (i.e., frequency = 3X10!°/X c.p.s.), 
A = [ A? 
(1.2) , = 1- 4a: 
1 = length of antenna in centimeters, 
Xo = distance of antenna from the central line of the wider face of 


the guide in centimeters, 
Zo = distance of antenna in front of plug in centimeters. 


The arrangement is shown in Fig. 1. 
§ £ 





CROSS SECTION OF GUIDE 


U 
5 
L 2 


LONGITUDINAL SECTION OF GUIDE 


Fic. 1. 
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As a particular result, we place the antenna symmetrically (xo = 0) and 
take / = \/4 and the following numerical values: 


(1.3) a = 5.5cm., b = 4.5 cm., A = 10cm.; 
then the resistance is 
(1.4) R = 185 sin® (0.262 zo) ohms. 


To match such a guide to a coaxial line with an impedance of 70 ohms we 
should take ae 
70 


(1.5) sin (0.262 zo) = iss = 0.615, 2 = 2.53 cm., 9.46cm.,... 


This gives the possible positions at which the plug should be set. 


2. The Lattice of Dipoles 
We use the method of images. This method is available only for the fol- 

lowing sections: 

(i) Rectangle, 

(ii) Equilateral triangle, 

(iii) Isosceles right-angled triangle, 

(iv) 30°—60°—90° triangle. 
The rectangular section is the simplest, and the only one of these sections 
actually employed; we shall confine our attention to it. 


The image of a dipole with respect to a perfectly conducting plane is well 
known. It is found by taking the geometrical image, and then reversing 
direction. The dipole and its image give a field which satisfies the condition 


(2.1) E, « ©, 


oa 
over the plane, E, being the tangential component of electric intensity. 


¥ 
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Let us first consider the case where the rectangular guide has no plug, so 
that it extends to infinity in both directions. Suppose that there is a dipole 
inside the guide. It is easy to find a doubly infinite system of images (all 
external to the guide) so that the lattice consisting of the dipole and its images 
gives a field satisfying the condition (2.1) over the walls of the guide. If we 
now add a plug, the condition (2.1) will be satisfied over it if we add to the 
above lattice its image in the plane of the plug. 

Let us take the origin at the center of the plug (Fig. 2), the axes of x and y 
parallel to the edges of the rectangular section and the axis of z along the axis 
of the guide, pointing towards the open end. Let a and 6 be the lengths of 
the edges of the section of the guide. 





Fie. 3. 


We place a dipole with components (Az, Ay, Az) at (x0, Yo, 20), these 
components being complex in general. Fig. 3 shows part of the lattice in 
the plane z = 29. The arrows are to be regarded as oblique to the plane of 
the paper and equally inclined to it. The signs indicate the components: thus 


(+, +, +) means (Az, As; A,), 
(+, -—, —) means (A.z, —A,, —A,). 


If we let w and » take all integer values, positive, zero, and negative, we 
may represent the components and positions of all the dipoles in the plane 
z = 2, as follows: 
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(4.;, a5, #p) at (xo + 2ua, yo + 2vb, 2), 
| (Az, —Ay, —Az) at (—%0 + (2u + 1)a, yo + 206, 20), 
De eid rs sill at: ame ies chemi 
(—Az, —A,, Az) at (—xXo+(2u+1)a, —yot(2v+1)), 20). 


The images with respect to the plane of the plug are given by this table on 
changing the signs of Az, A,, and 2. 

The lattice shown in (2.2) is a superposition of four basic lattices, in each of 
which all the dipoles are the same. The periods of all these basic lattices are 
(2a, 2b). Instead of calculating at once the field due to the whole system of 
dipoles, we first calculate the field due to the basic lattice shown in the first 
line of (2.2). The fields of the other basic lattices then follow by simple 
substitutions; adding all four fields, we get the field due to the system (2. 2). 
The field due to the plug is then easily obtained and added to the field already 
found. 

The transformation by which we pass from the field of one basic lattice 
to that of another is given by the following table, read downwards: 


Basic lattice Components Base-point 

I An. Ay, A, Xo, Yo, 20 

| II A,, —A,, —A; —xo +4, Yo, 20 
(2.3) 

| III —A,, Ay, —As Xo, —Yo + b, Zo 

| 

{ IV —A,, —A,, A; —xo ta, —yo t+ b, 20, 
The co-ordinates of a general dipole in the basic lattice I are 
(2.4) Xyy = Xo + 2ua, ¥,, = Yt 2vb, 2 = 2, 


where yp, v take all integer values, positive, zero, and negative. 
Let X be the wave length of the radiation (i.e., frequency = c/X), and let 
k = 2r/X. The time factor is e~**. We shall use an accent on a symbol 
containing the time factor, and drop the accent when the time factor is omitted. 
As usual, physical quantities are to be found by taking the real parts of the 
complex expressions. 
The Hertz vector at the point (x, y, z) due to a dipole A at (x’, y’, 3’) is 
= 
~ > | —> A eikr 
a) e = oe, $b =7-— 
where r is the distance between (x, y, 2) and x’, y’, 2’). Thus the Hertz 
vector due to the basic lattice I is 
oe 


— A @ S oetkruy 


(2.6) d = = 2 


Ar ure oO vr=—-@ r 
where r,, is the distance between (x, y, 2) and X,,, Vuy+ 2yy)- 


’ 





’ 
uy 
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3. The Field Due to a Basic Lattice: Discontinuous Method 


The idea used here is based on the familiar theory of a diffraction grating. 
It is known that a finite lattice composed of oscillating dipoles sends out 
pencils of radiation in a discrete number of directions. The field far from the 
lattice and very near to the center of a pencil of radiation can be represented 
by a plane wave; the field outside these pencils of radiation is practically 
nonexistent. Without going into mathematical details we can foresee what will 
happen if the lattice becomes larger and larger. All the discrete pencils of 
radiation become wider and wider, finally covering the whole space. The 
regions where the radiation was practically nonexistent shrink to nothing. 
The final radiation at a point, arbitrary but far enough from the lattice, 
can be obtained by taking into account all the plane waves traveling in all the 
admitted directions, given by the geometrical structure of the lattice. This 
is the principal idea which we shall work out here mathematically. 

Let us consider first 2 patch, that is a finite basic lattice, consisting of 
dipoles whose co-ordinates are given by (2.4), where 
(3.1) —-M<u< M, -N<vEN. 


Ss 

Let a, B, y be the direction cosines of the line going through the center of 
a radiation pencil and x, y, z the co-ordinates along this line. Then, at dis- 
tances great compared with the dimensions of the patch, we have, in the 
first approximation, 
(3.2) Ty = 1 — 2waa — 2vbp, 
r being the distance from the center of the patch; and the Hertz vector due 
to the patch is (to this approximation) 


-- “A em a en Sp tikeBs 
(3.3) OM, N) = 5 (3 3 )(\2, ). 


From (3.3) we obtain the possible direction cosines of the radiation pencils 
if Mand N are large, which we shall assume throughout. These are: 


(3.4) a= Am; B = A. Y¥ = Tans 


where 
5 — ae a oie a= 
G5) Geet. Bits te 1 _ a(S + a) , 


m and n being integers, positive, negative, or zero. Since y must be real, the 


condition for m, n is 
a. 4 
3.6 —+—< =: 
(3.6) a2 & «2 
the last inequality, together with (3.5), defining the directions of all the radia- 
tion pencils, or, in other words, all the (m, ) pencils. The Hertz vector along 


the center of the (m, 1) pencil is 


(3.7) _ ) = Ache 
o (M,N) ==" Mw. 
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This result is based on (3.2) which may be rewritten for the present choice 


of points of observation 


(3.8) Ymnyuv = Ymn — 2UAQ2m — 2VvbBn 
a 2uma _ 2vnt 5 
er k k 


The result as expressed in (3.7) is not very satisfactory since it contains 
M and N characterizing the dimensions of the patch. We can improve (3.7) 
by taking instead of (3.8) the rigorous formula 


(3.9) mn. = (Pan + 42a? + 40°b? — ArmnldQm — 4fmnvbB,)? . 


We see that (3.8) is the first approximation to (3.9), the second approximation 
being: 
(3.10) Tmna.uv = mn — a (um ae vn) 
2 m1 n>1? 2uvmnt? 
+ 2 (1 SE)+ (1-8) - 2] 


mn 





When this is substituted in (2.6), we get for the Hertz vector along the center 
of the pencil: 








a eitran a 
(3.11) Pmn (M, N) 4r mn be-M v=-N 
2ikT oo (4 _ mn? 8 ( - oR) - | 
exp [24 E ‘a (1 =) + yb? {1 be pk? ’ 


which is more rigorous than (3.7). 

Let us now increase M and N to infinity, that is, we spread the pencil over 
the whole space. But in this case the pencils overlap and air (M,N- ~) 
ceases to represent the Hertz vector inside the broadening pencil. It will 
represent only the (m, 1) contribution to the total 4 where 

oy mo ~~ — 
(3.12) ¢= 2% XS Pm (M, N- ~), 


m= -mo n= ~-no 


My and N, being the greatest integers satisfying (3.6). 
We turn now toward the calculation of 


= — 
(3.13) Dn _ Pn (M, N- o), 


by replacing the summation in (3.11) by an integration. This leads us to an 
integral of the type; 


(3.14) / i eich+dn'+22indtdn, 
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which can be brought by a rotation of the co-ordinate system into the form 


(3.15) ( eat) ( erin) = 4(c'd’)4 [za + o] = (c'd') ur. 


But 
(3.16) cd’ = d — g’. 


Comparison with (3.11) gives: 


(3.17) ¢ = 2ka (1 aE) d = 2k (1 — Fe 








_ _ 2mnr’* _ i le 
g ~~ kb ’ g — “wae n = ‘fF 
The result is: 
-: iA 
= ———__ eiktmn 
+ on * a 
where 
(3.19) mn = (X — Xo) Om + (Y¥ — Yo) Bn + (2 — 20) Yun» 


and (x, y, 2) represents an arbitrary point in space. The plane wave as 
represented by (3.18) has no physical meaning. To find the Hertz vector at 
a point, arbitrary but far enough from the basic lattice, we must perform the 
summation indicated in (3.12). Remembering that, according to (3.5) 
(3.20) an = —Aa; Bw os —Bn, 


we have, because of (3.18), as our final result: 


— 
o = tA {eiteom 22 = COS Ratm(x — Xo) et rmle-0) 


(3.21) as aa 


+ 2 3 = cos kB, (y a Yo) etk%on(z - 20) 
On 


n=1 


COS ROtm(x — Xo) cos RBn(y — yo)et*Van(2-%) 





my) we 
+42 2% 

m=1 n=1 mn 

representing the Hertz vector due to the basic lattice I. 
We can now find the Hertz vector of the lattices equivalent to that of the 
wave guide. We have to apply our formula (3.21) four times, changing the 


components of A and xo, Yo, according to the table (2.3). In this way, we 
obtain the Hertz vector of the lattice in the plane z = 29. The Hertz vector 
of the lattice in the plane z = —2 is obtained by changing zo into —zo , and 
A,, A, into —A,, —A, in this result. The final Hertz vector obtained in 
this way is a summation of the Hertz vector for eight properly chosen basic 
lattices. Denoting now by ¢,, ¢,, ¢, the components of this resulting Hertz 
vector we find because of (3.21) and (3.5): 
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2A. no : 
a y : : : o,.) pikYons 
= y+ — sin nx sin nXo sin(RY onZ0)e 

kab |y=1 Yon 





bz 


mo 


es A ; ; 
+2 5 DY — cos my cos my sin nx sin nxXo 


m=1 n=1 Ymn 


sin (RY mnZ0) ek" , 





CA Nee 1... F : oe 
oy = —) YS —sin my sin myo sin(RYmo20)e*™ 











a kab aed Y mo 
(3.22) ‘ 
mo me 1 “ ‘ 
+23 DY —sin mp sin my cos nx cos nXo 
m=1 n=1 mn 
sin (RY mnZ0)et* , 
4iA, mm : 3 : ‘ 
¢: = oe sin my sin mo sin nx sin nXo 
kab m=1 n=1 Ymn 
cos (RY mnZ0) ek" 
where 
i . ¢ 
1x T TX 7 
j y ee ee ae le Yo ge 
| a 2 a 2 
(3.23) 
_ my 7 _ Vo v 
x b 2’ Xo b 2 


Equations (3.22) represent the Hertz vector for great z inside the wave guide, 
i.e., for 
—ta gx < ja, —3b Sy < hd. 


4. The Field Due to a Basic Lattice: Continuous Method 


Let us consider the basic lattice I in the plane z = 20, as in (2.3). Instead 
of discrete dipoles, let us take a continuous distribution over the plane with 
density o,so that a dipole o dx’dy’ is associated with the element dx’dy’, 
where x’, y’ are co-ordinates in the plane z = 2. Let p be a large positive 
integer; consider the expression 


(4.1) Cv = 422 cos””.X’ cos?” Y’, 


ab 


where K(p) is a constant and 
(4.2) X’ = 437(x’ — Xo)/a, Y’ = 32r(y’ — yo)/d. 


It is evident that this density will be very small except in the immediate neigh- 
borhood of the positions of the dipoles in the discrete lattice; the larger we 
make p, the more closely does this continuous distribution approach the dis- 
crete distribution. We choose the constant K(p) so that the total dipole 
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strength inside the period rectangle (2a, 2b) is the same in the continuous and 
discrete distributions, i.e., so that 


(4.3) [ [eevay mA, 


the ranges of the integral being 
XH -aSxwSxta w-b< W< w+ 5. 


Since 
xXot+a ) jr 
| cos??X’ dx’ = = | cos’??X’ dX" 
xo-a -4}r 
= (2p)! 
~ 74 Qe(pip 
we take 
“(py x 2rCb!)* 


By Stirling’s formula this is approximately 7p/4, but it is actually more con- 
venient to retain the exact value (4.4). 
The Hertz vector for the basic lattice, smoothed out as described above, is 


: = 1 eae, aes 
(4.5) @ = is oe" /r .dx'dy’, 


where 7 is the distance from the point (x’, y’, 29) to the point of observation 
(x,y, s). Let us drop a perpendicular from (x, y, z) on the plane z = zo, and 
take polar co-ordinates (R, 9) relative to the foot of this perpendicular (Fig. 4). 
Then 


(4.6) x’ = x+ Roos 8, y = y+ Rsin #9, 
and (4.5) gives 
s > AK), ao 
(4.7) @ =- al cwrsrar | cos??X’ cos” Y'd@. 
0 0 
On substituting ae 2» 
nS y’ 2imX’ 
cos?X’ = 2 ,~, % a m)é , 


and a similar expression for cos ?”? Y’, we have 


(4.8) | cos??X’ cos?” Y’ d6 
0 


1 p p 2 2 ox aR mcos@ — nsin@ 
= | ae eS a p p e2i(mX + nY) — (most + 6 J dé 
2M nSepccue tt + MIO +s : 
0 


where X, Y are given by (4.2) on removing all accents. Since for any real 
constants k,, k:, 


(4.9) l eitkicosd + kin) dQ = 2rJo(+/ki + R3), 
/0 
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where Jo is the Bessel function of order zero, (4.7) reduces to 


—- Be as 


Zz 


POINT OF OBSERVATION 





Fic. 4. 

where 
(4.11) Tos -/ 6M rR (Cacia, 

0 
(4.12) Cnn = 1/m?2/a? + n?/b?. 
Now 

R? = r — (g — 20), 

and so 
(4.13) | -/ e*T6(CnnV P a (2 = Z0)*)dr. 

|z-20] 


This integral is evaluated in Appendix A; we find 





(4.14) a 


eth Y¥mn\ 3-20] 3 
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where Ymx is defined by 


( RY mn = | VR = 3, | if Crn < k, 


(4.15) iced 
| RY mn = i| Ve — mn if Cun > k . 


When we substitute this in (4.10) we get for the Hertz vector of the smoothed 
basic lattice I 
> iAK() & & ( ® \( » \1 
= s > 
(4.16) ? 2 + kab ms-» n= -p ks ve ad (, 7 n) Ymn 


eiklam(x-x0) + Bn(y-9) + Ymn|2- 201] : 








where 
kan = mt/a, kB, = nw/b. 


If we go a long way down the guide, so that z — ©, only those terms in 
(4.16) survive for which Ymn is real, i.e., 


(4.17) Cun < b, 
or equivalently 

m? n? + 
(4.18) ae ee te 


This condition is the same as (3.6), and is in fact the cut-off condition given by 
the eigenvalue method for an empty guide. The condition limits the values 
of m and n to some finite ranges 


—m <m<m, —mo Ln Km. 
. . — . rT 
Thus in calculating @ from (4.16) we have only a finite number of terms. We 


now proceed to the limit p— © in order to pass from the smoothed basic 
lattice to the discrete basic lattice. Since, for fixed values of m and n, 


a 
om x be (, - J 
se i 8 ae ee tea 
prot (p+m)l(p—m)ipt+nji(p—n)l 4 
we get for the Hertz vector of the discrete basic lattice I a long way down the 
guide 





> 
<==) 1A mo no 1 4 
4.20) = y y eiklom(x - x0) + Bn(y- yo) + Ymn(2 -20)] 
(4.20) P = Shab mmm nam Yuu 


where 


a N/m? a 
(2) yee ft - Bt - FG +9). 


Equation (4.20) is the same as (3.21), and so the two methods lead to the 
same expression for the field at a great distance down the guide. From this 
point the argument proceeds in common. 


However, before finally leaving the smoothed lattice, let us write down 
the Hertz vector at any point inside the guide. To do this, we have to super- 
impose on (4.16) the corresponding expressions for the other basic lattices in 
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(2.3), and then superimpose on the result the expression corresponding to 

the images with respect to the plane of the plug. A simple calculation gives 
“at *K(p) & 2p 2p 7 

4.22 = — : at Or > ( eik( mx + Bny) 

( ) ? 24? + kab cay n i p +m p +n Ymn 


— ; ~ hs ; , 
[ Fumnye!* Ymn!2 - 20| + Gomnyetk Ymnl2 + 01] : 





where 

( Frans = —Gounys = Azle Hoe + (—1)metm™] [e-thBae — (— 1) sehAarn] , 
(4.23) } Finny = —Gimny = Aylenitamt — (—1)metkam] [e-ikBnso 4 (—1) eit Bnrr) , 
A, [e~tkamts — (— 1) meikam*o] [e—ikBnyo — (—1)mei®Bnye] , 


\ F ante = Giawds = 


5. The Hy Wave 


Let us suppose that the dimensions of the guide and the frequency are such 
that 


(5.1) 1 4 1 1 

yo ES + b 
These conditions are satisfied for example by the numerical values given in 
(1.3). Then the condition (3.6) or (4.18) is satisfied only by 


(5.2) m = —1, 0, 1, n = 0. 
Putting mp = 1, mo = 0 in (3.22) and writing 
: ? 
(5.3) YS ae = 0. 1 = 4a?’ 
we get for the Hertz vector a long way down the guide 
dz = @: = 0, 
(5.4) 2A Tx Tx 
} Lay Xo . A Jikyz 
= - cos —— sin kYZo cos es 
Pu kyab a ™ a 
The corresponding field is 
. 0” @ ¥ 
E =- —_—- = 0, 
" Ox dy 
‘ od, , s 
Ey oe k*dy 
oy- 
2kA , TXo . TX ikys 
= — cos — sin kyZ) cos—e ', 
yao a a 
(5.5) (> 
Eve = Oe 0, 
Oz0Y 
a C8 QkA WXo . TWX ikys 
H, = -4k Py = — —— cos — sin kyz0 cos — e ~ 
Oz ab d a 
H, = 0 
.. Od, 217A, TXo . - WL abr 
H, = —itk- Py = — cos — sin kyzo sin — e"”’. 
Ox ya*b a a 
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The time-mean of the flux of energy across a section of the guide is 


(5.6) F = real part of | | (E.H, — E,H,) dx dy 


RcA,Ay ety... 
oS cos* —— sin? kyzo. 
~ yab~ a 
The transmitted wave is the Hj) wave, as indeed we know otherwise by the 
eigenvalue method. It is the quantitative character of the present result 
that isnew. We note that: 
(i) The radiation depends only on the component A, of the dipole, 
i.e., on the component parallel to the shorter edge of the cross 
section. 
(ii) The radiation depends only on two (xo, 20) of the co-ordinates of 
the dipole. 
a 
Acurrent I’ = I exp (—ikct) flowing in an element dso (dxo , dyo , d2o) of wire 
with direction cosines L, M, N is equivalent to a dipole with components 
ilL iIM iN 
2 dso, Ay = <- Gla, 44 = — dso. 
‘If the current element belongs to an antenna in a rectangular wave guide 
with a plug, the field a long way down the guide is given by (5.5) when we 
substitute from (5.7). 
If the antenna is very short (length 7) compared with the wave length, and 
the current J is regarded as constant along it the output is by (5.6) 
Mere TXo 


(5.8) F = cyab ‘ -os" — sin? kvyzo. 


) A, = 


~I 


(S. 





Since the time-average of the square of the current is }//, the radiation 


resistance is 

a F 2M??? B60. «- 

(35.9) R = — = — cos? — sin’kyZo 
17] cyab a 


in Heaviside units. Since 


(5.10) 1 Heaviside unit of resistance = 47 X 9 X 10" ohms, 
the radiation resistance of the short antenna is 
z MIP i 
(5.11) R = 2407 — cos? — sin? kyzo ohms, 
yab a 


the lengths being measured in centimeters. 
For a finite antenna the surviving parts of the field are by (5.5) 





_ a a , me ies 
i cyab “ _" 
5 _ aT TX skys 
2 a mx 
(5.12) i, ~~ ee e's 
2 J e x 1RVZ 
H,.= = Ps A OR . el 
ckya’ b-~—Ca 
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where 

(5.13) J= / I cos = sin ky20 dyo, 

the integral being taken along the antenna. Accordingly the output is 
; JJ 

(5.14) a 


If Jy is the current entering the antenna, the radiation resistance is 


F QJ 








(5.15) ice =— = ——— Heaviside units 
41 lo cyablolo 
= 200r -L- IJ ohms. 
yab Iolo 


When this is calculated in any particular case the current Jo cancels out, and 
the result depends only on the shape and position of the antenna and the point 
at which it is fed. 


6. Straight Antenna Parallel to the Shorter Side of 
the Cross Section 


When the condition (5.1) is satisfied, the only radiation transmitted by 
the guide is the Hijo wave. Any form of antenna will transmit this wave, 
except an antenna lying entirely in a plane perpendicular to the y-axis. But 
it is only the component of current parallel to the y-axis that contributes to the 
radiation, and the simplest plan is to use a straight antenna parallel to the 
y-axis. 

Let us consider the arrangement shown in Fig. 1, where the length of the 
antenna is / and it is fed at the point where it enters the guide. Let J be the 
current entering; then according to the usual rule of sine-distribution the 
current along the antenna is 


sin k(l — $b — yo) 


(6.1) cee sin Bl 
Then (5.13) gives ‘ 
pl-3b | R(I 1p ) 
si TXo sin — 20 — Je 
(6.2) J =: Jg'cos : sin bre | aa dyo 


= I cos oe sin kyzo : tan $l, 


and so by (5.15) 


= 1 ge TOs 21 
(6.3) R = 240 ™ Byab cos’ ~~ sin kyzo tan? $k/ ohms 


2 
* x cos? = sin? ouyse tan? xt ohms. 
w yab a nN 
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Let us now place the antenna centrally, so that x» = 0, and at the same time 
put / = \/4; we get 


(6.4) _ 60 


~ a yab 


* 27 Zo 
sin? va ohms 








Finally we take the numerical values given in (1.3) and obtain (1.4); the 
matching condition (1.5) follows at once. 


It is clear from (6.3) that there are three different ways of matching the 
radiation resistance of the antenna to the impedance of a given line: 


(i) By adjusting the position of the antenna relative to the central line, 
i.e., by changing xo; 
(ii) By adjusting the position of the plug, i.e., by changing 2 ; 
(iii) By adjusting the length of the antenna, i.e., by changing /. 
Until such time as reactance has been calculated, it is not possible to decide 
which of these methods is the best. 
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APPENDIX A 


We shall now evaluate the integral occurring in (4.13), namely 
(A1) I(f) = [ e#® Jy (C VP — &) dr, 
£ 


where ¢, k, C are real and positive. 


Differentiation with respect to [ gives 


ee ae TCV = $) 
(A2) dt e ct | e Joe dr. 


We integrate by parts, writing Vr? — [ = R; 


co , @ ikr 
(A3) C | gir LAER) oy in — Ji(CR)d(CR) 
§ r= 


t c 


eikr ? ' d fer 
a E rcr) |. a | . ) Jo(CR)dr. 


The first term here is —e*5/f. Hence (A2) becomes 


es d etkr 
. I 5(<) Jo(CR)dr. 


(A4) 


“rl = 
S/S. 
eh mm 
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Comparing this with (A1), we observe that the application of the operator 
s 
¢ df 
to the iniegral has the effect of inserting the operator 
ai 
drr 
in front of the first part of the integrand. The fact that that part is actually 
e*" is not involved in the argument, except in so far as it secures convergence. 
Thus, applying the same operator to (A4), we get 


(AS) : f. (3 a7) ss ie {3 £ (=) Jo(CR)dr. 


Let us now go back to (A1) and put, by virtue of the Bessel differential 
equation, 


(A6) Jo(CR) = —Ji’(CR) — Jo(CR)/CR. 
Thus 
(A7) Ig) = - of e"(R/r)Jo'(CR)d(CR) — cf e*"Jy(CR)/R.dr. 
r=¢ c 
On integrating the first integral by parts, we obtain 
a , d aie A 
(A8) (gf) = Cf (5. (c# *) _ I J((CR)dr 


Js 
= C1 i on(@ “ 3 )RIGCR)dr 
c 


és ony eu(S = 3 )RJW(CR)A(CR), 
rat 


r 


and further integration by parts gives 

(A) CUS) = i s, {en(- B+ 4) @- )} JCR)ar. 
; dr prey 

This integral splits into two parts. The first is 


ihe 2 eis 
(A10) I = ie (- ik + )} Jo(CR)dr 


ro @ ikr 
| e*rJ\(CR)dr +f 5 (F)s(cRdar 
c £ 


1d 
d 


Il 


~ 


’ 


RI + 


yi 
“~~ 
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by (A1) and (A4). The second part of the integral is 


“ of4 en(- 3 s+ 1) solCR)dr 
= fz a {3 £(¢ ")} Jo(CR)ar 


=H? Re) 


by (A5). Thus (A9) becomes 


(A11) 


1 dI d/1dI 

A12 CE = Be SS He a) 
oe Witeagt! HER 
or 
(A13) 7 

oa ee (2 — C*?)I = 
Hence 
(A14) 1(f) = AetVP-4 4 Be-tV#-2, 


where A and B are independent of ¢. For definiteness, and without loss of 
generality, we interpret the radical in the following sense: 


VP —-C = |\VP—-C| if k>C, 
VE —C =il\VPF—-C| if k<C. 


Now putting ¢ = 0 in (A1), we have a well known formula (Gray, Mathews; 
and Macrobert, Bessel Functions, p. 65) 


(A15) 


(A16) Iz20 = f e*" Jo (Cr) dr 
0 
‘ii ae 
VE-CG 


the radical being interpreted as in (A15). Also by (A4) 


(A17) (=), = lime {3 (=) ICV Pr — €)dr 


dé f—-0 


. mf SE =) Jo(Cr)dr 


= lim, ¢ { - = Jo(Cf) - cf  J5(Chdr] . 
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Since the last integral has a finite limit for £[ — 0, 


we get 


— (zt), ae 


The values (A16), (A18) enable us to determine A and B in (A14). We 
obtain 


(A19) I(f) = | e* J (CVP — @)dr = a eitVERe 





VRP - C 


the required result; the interpretation of the radicals on the right hand side 
is as in (A15), which agrees with (4.15). 

We may remark that our result may be obtained formally by changing 
reals into imaginaries in a known formula; cf. Riemann—Weber, Differential- 
gleichungen der Physik, vol. 2(1927), p. 546. Such a procedure is not, of 
course, mathematically valid. We have not been able to find our result in 
the literature of Bessel functions. 


APPENDIX B 


The following treatment of the problem considered in this paper does not 
make use of the method of images, and offers interesting possibilities of 
generalization. In this, we find the steady-state solution corresponding to 
an assigned volume distribution of current throughout the wave guide, and 
deduce the field due to a current element (dipole) as a limiting case. 


We use the same notation as before, except that it is now rather more con- 
venient to take the origin at a corner of the plug-face than at its center. In 
Fig. 2, the origin must therefore be supposed transferred to the bottom left- 
hand corner, so that the equations of the faces of the guide are x = 0, x = a, 
y¥=0,y=b. 

aa 


Maxwell’s equations for a volume current density 7 may be written, in 


Heaviside units, 
ao — — — cack 
(B1) curl E = tkH, curl H = — ikE * - 


We shall, for the purposes of calculation, suppose that the wave guide 
extends to infinity in both directions along the z-axis, and eliminate the plug. 
We then assume that j,, j, are odd functions of 2, and that 7, is an even 
function of z. It is then evident that the boundary condition (2.1) will be 
satisfied over the plug-face (as may be verified from the solution to be 
obtained). 


We now make the following assumptions for the field components, which 
obviously satisfy the boundary conditions (2.1) on the faces of the guide: 








CHIEN ET AL.: CONTRIBUTIONS TO THE THEORY OF WAVE GUIDES 89 
E, = 2 ASn (z) cos (Cnx) sin (D,y), 
E, = 2 AS. (z) sin (Cmx) cos (Dny), 


|= 2 AS, (2) sin (Cux) sin (D,y), 


(B2) 


RB: « = Bon (2) sin (Cnx) cos (Dny), 
H, = 2 B®. (z) cos (C,x) sin (D,y), 


H, = % BE) (2) cos (Cux) cos (Dry), 
where C,, = mm/a, D, = na/b, the summations running from zero to infinity. 
—_ 
For the components of 7 we now assume 


je = 2 582 (@) cos (Cux) sin (Day), 


(B3) iv = & js (2) sin (Cnx) cos (Dny), 
je = Z 5 (@) sin (Cux) sin (Dyy). 


Substituting (B2) and (B3) in the Maxwell equations (B1), we have six linear 
equations of the first order to determine the six functions A%» (sz), 
B® (z) in terms of the functions 79), j,, 72. The solution of these equations 
is facilitated by noting that from (B1) follows 


=> 1 —> 
(B4) (V27+h) H = — . curl 7. 
Taking say, the x-component of this equation, and using (B2), (B3), we find 





2R() (2) 
(BS) ot t+ BR = + (4 - im), 
where 
(B6) Unn = k? — C2 — D2. 


For definiteness, we suppose the positive root of (B6) is taken for tmn if u2,, > 0 
and the positive-imaginary root if ui, < 0. The solution of (B5) which 
represents outgoing waves, or damped waves, for large positive or negative 2 
is readily found, by the method of variation of constants, to be 


(B7) BQ = _ [ 00 +5 |? fo» ) (z)dz, 


(3) 








where QB = e “ar (ites JQ — D jG 


and the bar denotes, here and in what follows, the complex conjugate (the 
functions j“;2,2 being treated as real). This gives outgoing waves at infinity 
if Umn is real, and damped waves if umn is positive imaginary. 














90 CANADIAN JOURNAL OF RESEARCH. VOL. 27, SEC. A. 


In a similar manner, the y- and z-components of (B4) give us B®, and B®), 
and the second of (B1) then gives us A%,, A®,, AS. We thus find: 


§ Z 
eiUmn2 
Bas = 5 | (2)dz + 5 
-2 





a f'2Q(2)dz, (a = 2,3), 











27CUmn 2icum 
; 2 : a 
e'Umn? e7 mn? be 
(B8) A® = x | (8) oe ® ()dz — + ig Sea 
as i 
(B = 1,2,3; 6g3 = 0 unless B = 2. 633 = 1), 
where 
= emit (Ca JQh — ittmn J), 
= etme (D,IR — Cuih), 


gi = entume (C2 — R)G2 + CaDrji2 — iCntimaG 2), 
g2 = em itm ice Dj ote (D2 - Rk?) 72, i 1Dn Nun Jomls 
& = etn [$Cnttan IQ + aD. nttnn Jn + (un — a) Inn : 


In deriving (B7), (B8), it has been assumed that r, and hence the functions 
jo,2, tends to zero at infinity sufficiently rapidly, and also that 7, , j, , and 
hence j{},”, vanish for z = 0, so that these functions are continuous for z = 0 
If this latter condition is not satisfied, other terms appear; but we can evidently 
suppose it satisfied, since we are going to proceed to the limiting case of a 
current element. 


(B2), together with (B7) and (B8), now gives a solution of the equations 
(B1), as may be verified by direct substitution, and all conditions are satisfied. 


Suppose now that the guide transmits only the Hi) wave; then from (B6) 
we see that all the wm, are pure imaginary except %9. From what has been 
said above, it then follows that all the functions AQ, ..., BS, represent 
damped waves except those with m = 1,n = 0. Hence, from (B2), (B7), (B8), 
and since the j& — 0 as z > ~, we find for the field, for z large and positive 


@ 
—-k_, : : -. We 
E, = emo | 7{2e-mo? dz. sin —» 
2cu10 a 


co 


| 





(B9) he aa “Ey, 


@ 
Wt. . Wx 
H, = emo? | 7{De- iw dz. cos — » 
2acuyo a 


-@ 





This represents an Hi) wave. From (B3) we have 


a b 
2 ; ? 
(B10) iB) == dx i dy. x(x, y, 2) sin. 
0 0 
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To obtain the field due to a current element of components (ILdso , JMdsp , 


pat 
INdso) situated at the point (xo, yo, 20), we now suppose that 7 tends to zero 
everywhere except in the immediate neighborhood of the points (xo , yo , 30) and 


(xo, ¥o, —20). At these points 7 is to tend to infinity in such a way that the 
total components of current in any small region surrounding (xo, yo, 20) are 
(ILdsy , IMdsy), INdso), and the total components of current in any small 
region surrounding (xo, yo, —2Z0) are (— ILds), —IMdso, +INdso). The 
introduction of the singularity at (xo , yo , —2») is, of course, necessary in order 


—s 
to preserve the odd or even character of the components of 7. This case can 
be represented symbolically in terms of 6-functions by putting, in our final 
results, 


je = ILdso6(x — xo) 6(y — yo) [6(2 — 2) — d(z + 20)], 
(B11) jy = IMdso6(x — x0)d(y — yo) [6(2 — 20) — 5(z + 20)], 
je = INdso6(x — x0)5(y — yo) [6(2 — 20) + O(2 + 25)]. 
Using (Bio), we then find, for this case, 
(B12) [iv e~“st dg = =* sin = sin (u9%0) 1M dso. 


Inserting this in (B9), and putting u:.0 = yk, where y is as defined in (5.3), 
we finally have for the field at a large distance down the guide 


( _ 2tIMdso . 4x0. eo 

| EB, = ~yabe sin — sin Vk20 sin an e '. 
(B13) | tw we 

| 

| 


2mIMdsy . xo . + Ee 
——-_—_— sin — sin ykzo sin = etrke 


H, = 
( 


Comparing with (5.5) and (5.7) and allowing for the change of origin, we 
see that this field is identical with that given by the method of images. The 
rest of the calculation now proceeds as before. 


yka*be 


It may be remarked that use of (B11) in the general formulae (B2), (B7), 
(B8) gives the field at amy point in the guide, expressed as a double Fourier 
series (except in the plane z = 29 , where the representation obviously fails). 


PART II. A GENERAL METHOD FOR CALCULATING THE IMPEDANCE 
OF AN ANTENNA IN A WAVE GUIDE OF ARBITRARY CROSS SECTION 
By A. F. STEVENSON 


1. Introduction 


In this part, it is shown that the problem of finding the radiation resistance 
and reactance of an antenna of arbitrary shape in a wave guide of arbitrary 
cross section can be solved if the current distribution in the antenna is known, 
and if the eigenfunctions which determine the free propagation in the guide 
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can be found. In the case where only a single E- or H-wave is transmitted, 
a simple closed expression for the resistance results. The reactance, on the 
other hand, requires the summation of an infinite series, the calculation of 
which would be somewhat tedious, though quite feasible, in any given case. 
The reactance is discussed in Section 4. Incidentally, the field at any point 
in the guide is obtained in the course of the work. 


The method used is a generalization of one of the three methods previously 
given for calculating the field due to a dipole in a rectangular wave guide in 
Part I. Its great advantage, as compared with previous special methods for 
rectangular and circular guides, is its generality and comparative simplicity. 


In a recent report, Chu (2) has given formulae for the resistance of antennae 
in rectangular and circular guides in particular cases, and has stated that the 
method used is applicable to guides of any cross section. He gives no 
indication of the method, however, beyond stating that it involves the use 
of orthogonal functions. Actually, some of Chu’s formulae are not correct. 


2. Calculation of the Field 


We consider a cyclindrical wave guide closed by a plug, both guide and 
plug being perfectly conducting. We take origin in the plane of the plug-face, 
with the z-axis along the length of the guide (Fig. 1). We use Heaviside 


x 






s 
ZA Ty YoZ) 
XoYoZo) 


: 
Fic. 1. Wave guide with antenna of arbitrary shape and current distribution. 


units, and suppose a steady state in which all quantities vary with the time 
according to the factor e~'!, this factor being omitted (the electric field is 


— 
thus the real part of Ee~*, etc.). 


We first recall the solutions which give the free propagation in a guide 
extending to infinity in both directions. For E-waves, we have to solve the 
eigenvalue problem 

ey , &y ig otk 
(1) ani tay t HY = " inside S, 
y = 0onC, 
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where S denotes the interior of the cross section of the guide, and C denotes 
the boundary of S (Fig. 2). Let the eigenvalues of this problem, arranged 


b 


re) x 
Fic. 2. Cross section of guide. 


in ascending order of magnitude, and the corresponding orthogonal eigen- 


functions, be i i i 


yY = hh, r,... 


We suppose the eigenfunctions normalized so that 


[ eas = 1. 
s 


The case of degeneracy can be included by supposing that some of the y’s 
may be equal. The field for an E,-wave, propagated (or damped) in the 
positive z-direction, is then given, apart from an arbitrary constant factor, by 


| B= KOBE. By = KOE, By = wads. 





(2) 





be e- an ad" aWn a 
| Hy = — ikfs(2) Git Hy = ifs) Got» He = 0, 
where 
k = w/e = 2x/r, f(s) = &, 
and 
(2') un = VE= A. 


In (2’), the positive square-root is to be taken if u, is real (propagation), and 
the positive-imaginary if u, is pure imaginary (damping). 
For the H-waves, we have to solve the eigenvalue problem 
OV , SY 4 uw = 0 inside S, 
0x" Oy 
av 
Ov 
where v denotes the outward normal to C (Fig. 2). Let the eigenvalues, 
arranged in ascending order of magnitude, and corresponding eigenfunctions, be 
M = M,, Me,..... ; 
Vv = YV, ’ WV, ee De ae ak 





(3) 
= Qon C, 
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As before, we suppose the eigenfunctions normalized so that 


[viecas = 1, 
s 


and degeneracy can be included as before. Then the field for an H-wave is 
given by 














E, = E, = -% E, = 0 
(4) } av, , 
H, = F,(z) —* H, os MAF. , 
where 
F,(z) = etU=, 
and 
(5) = VF — MM, 


the same convention as to aa sign of the square root being used in (5) asin (2’). 


We now consider the problem of finding the field due to a current element 
in the wave guide. For this purpose, we suppose the guide extends to infinity 
in both directions, and eliminate the plug. We then suppose that there is an 


—_ 
assigned volume distribution of current, of density 7, throughout the guide 
and find the field due to it. We then find the field due to a current element by 


—s 
proceeding to the limiting case where 7 vanishes except at the element and 
its electrical image in the plug-face. 

The Maxwell equations for the case considered are 


aia ae ao ale —_ ~ o- 
(6) tk H = curl EZ, thE = —curl H += j. 


We suppose that 7 is continuous and vanishes at infinity sufficiently rapidly. 
Bearing (2) and (4) in mind, we now assume for the field a superposition of 
E- and H-waves of the form 


E, = z[ a, (2) oy + ik A,(2) | 

















dy 
ei Z| an(s) a — ik A,(2) am 
7) E, = & biz) Yn, 
i, = z| -it ca(2) aa + C,(2) | 
H, = | it cals) aye + Ci = |. 





H, = 2 D,(z) Vv, ’ 


where the summations run from 1 to ©, and where a,(z) , ---, D,(z) are six 
functions to be determined. The field (7) satisfies the boundary conditions 
of the problem, and it can be proved that the most general such field can be 
expressed in this way. 
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— 
For the components of 7 , we assume eae 


* js se z[ pa) 3 (a = =| 
1 n 
[ pate) 9 (z) i) , 
rn(2) ast ; 
We now note that by differentiating the first two of (8) with respect to x and 


(8) “i= 2 
y, and combining, and using the differential equations (1), (3) satisfied by the 
functions y,, V, , we can replace (8) by 








are 
i" 


aM 


wa 
| 2 i, = — Bat Pale) Vos 
(8’) ! * fp = — ik EME quls) Va, 
. 
> ye ~ rn(2) Wns 


where 
j= e+ Se, jp = Se - 
x oy oy Ox 
We shall refer to (8’) as the transform of (8). In a similar way (by differentia- 
tion and combination of the x- and y-component equations), we construct the 
transforms of (7) and (6). 

We now substitute the transforms of (7) and (8) in the transform of (6). 
We then obtain six equations, in each of which only one of the sets of functions 
Vi, ve, --- or Vi, Ve, --- occurs. We can then equate coefficients of each 
y, or WV, on both sides of the equations. We obtain in this manner three 
differential equations of the first order for the functions a, , 5, , c, and three 
for the functions A, , C, , D, , the functions p, , gn , 7, being regarded as known. 
These differential equations are quite simple and are easily solved by the 
method of variation of constants, the arbitrary constants in the solution being 
determined so that only outgoing, or damped, waves exist at z = +o. 
We thus obtain 


an(2) = ch(z) — 5 Pal), bn(z) = p2 cr(z) — Or 


eitns - at en itnz 
Cn(2) 2 en | (iUnPn — r,)dz — a a x 
(9) / "= (tunPn + 7n)d2 , 
1 
A,(2) = ie D,(2), C,(z) = 2 Di (z) , 





2 s fz o 
D(z) = — ee eve | eis o, ds — oe ive f Ut g, dz. 
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The constants u, , U, are those defined in (2’), (5) with the same convention 
as to the sign of the square root. 


— 

The functions p,, gn, 7m are found in terms of the components of j by 
applying to (8’) the usual process for determining the coefficients in expansions 
in orthogonal normalized functions. We thus have the field due to any 


a 
continuous distribution of current 7. We pass to the case of a current element 


on 
of strength Jdsy situated at the point (xo , Yo , 30), in the presence of the plug, by 
putting, in accordance with the principle outlined previously, 

jz = I, dso 6(x — xo) b(y — yo) [6(2 — 20) — (2 + 20)], 
(10) jv = I, dso 5(x — x0) 6(y — yo) [6(2 — 20) — O(2 + 20)], 

jz = I, dso 6(x — x0) O(y — yo) [6(2 — 20) + O(s + 20)], 
the notation of ‘‘delta-functions” being used. Finding first the functions 
Pn» Qn, % by means of (8’), substituting in (9), and then using (10), we thus 
find, forz > Oandz ¥ 2, 
1Un 4 


ll 





a,(z) = Pane P,(X0, Yo)filtln » 2, 20) — ke Ri(Xo , Yo)fo(un , 2, Zo), 
: 2 
b,(z) = — c Pi(xo, Yo)fa(tn , 2, 20) — ew, Ra(tos yo)falte » 8, 80), 


1 
n\Z =—), ’ 
cx(s) = 35 dale) 


(11) : 
A, (3) ~ cU,M.2 Qn(Xo . yo)fi(Un » 2, 20), 
1 
C,(2) = cM.2 Qn(Xo ’ yo)fs(Un » %, Zo), 
Ds) = M2A,(¢), 
where 
Oi IWn(Xo, Yo) dWn(X0 , ve] 
Palxo ’ ¥o) = [7 Oxo —s + ir . dso ’ 
(12) ; = [ OV (Xo , Yo) = OV, (xo, v0) 
Qn(Xo, Yo) = | 1, dy I, in dso, 
R,(xXo ’ Yo) = L.Wn(xo ’ yo)dso ’ 
and 
Si(u, 3,20) = ef sin uz, 3% > 20, 
= é*% sin uz,3 < 20, 
fo(u, 2,20) = e'“* cos u2%o,2 > 2, 
= te'“% sin uz, 2 < 20, 
(13) fs(u, 2,30) = te’ sin uzo,2 > 20, 
= e*% Cos uz, 2 < 20, 





e*2 cos uz, 3 > 20, 
= e*0cos uz,3 < 2%. 


Sa(u, Zz, Z0) 
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Equations (7), (11), (12), (13) now give the field at any point in the guide 


due to a current element T dsp at (Xo, Yo, 20) in the presence of the plug. It 
will be observed that the functions fe , fs are discontinuous at z = z). This 
is due to the nature of the field in the immediate neighborhood of the element; 
in the plane z = 20, the series (7) are not valid in general. But this is of no 
practical consequence. 

To find the field due to a thin antenna with a given distribution of current of 
linear density T(x , Vo, 20) (Fig. 1), we have merely to write I = T(xo » Vo. Zo) 
and integrate the expressions for a,(z), ---, D,(z) given by (11), (12), (13) with 
respect to so along the antenna, ds») being an element of length. We can draw 
the following two general conclusions: 






(x= Xo »¥=Yo ) 


Fic. 3. Longitudinal antenna. 


(1) A longitudinal antenna (Fig. 3) can only generate E-waves, whereas a 
transverse antenna (Fig. 4) can generate both E- and H-waves. This 
fact is already doubtless well known to experimental workers. 


‘J ‘Kz-20) 


Fic. 4. Transverse antenna. 


(2) When a transverse antenna is used for generating an E,,-wave, the effect 
is a maximum, other things being equal, when the antenna coincides 
with one of the orthogonal trajectories of the curves Y, = constant. 
The effect is zero when the antenna coincides with one of the curves W, 
= constant. 

(3) When a transverse antenna is used for generating an H,-wave, the effect 
is a maximum when the antenna coincides with one of the curves V, = 
constant. The effect is zero when the antenna coincides with one of the 
orthogonal trajectories of the curves VY, = constant. 
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3. The Formulae for Radiation Resistance 

To find the resistance, we need merely calculate the mean flux of energy 
at z = + by Poynting’s theorem, and equate this to the mean rate of working 
of the applied voltage. The method is sufficiently well known, and we omit 
the details of calculation. In considering the field at infinity, only a few terms 
(usually only one) of the infinite series in (7) survive, all the other terms being 
damped out. We consider the following two cases: 

(1) The dimensions of the guide (and manner of excitation) are such that 

only the E,-wave is generated. We find 


if (« oy. +I, oe) sin (u,z) ds 


nck / T.Wn COS (Uns) ds 


An 














oes 
| Zo|? 
where ¥, = wW,(x, y), Jo is the current at the feeding point, and 


ku, u,/R 


ft 1 (u,/kP 
We note that u, = 27/X’, where ’ is the wave length in the guide, so that 
u,/k =X/X'. Theintegral in (14) is extended over the length of the antenna 
(we have dropped the suffix 0). We recall that y, is a normalized eigenfunction 
(for definition of un, Wn see beginning of Section 2). If u, is a degenerate 
eigenvalue of (1), (14) must be replaced by the sum of such expressions for 
each independent and orthogonal eigenfunction y, . 

(2) The dimensions of the guide (and manner of excitation) are such that 

only the H,-wave is generated. We find 


19 
1 av, eS ie = 
if (- ay -—TI, =) sin (U,z) ds To|? ’ 


casi a 
"= U.Mi ~ 1— (U,/k)? 


n 








(15) (R)ohms = 12078, 








where 


The same remarks apply to the case of degeneracy as in (14). (For definition 
of M,, Un, Vn see beginning of Section 2.) 

We conclude that for a longitudinal antenna lying along the line x = xo, 
y = yo, and generating an E,-wave, R varies as [Wa(xo, yo)]?; while for a 
transverse antenna lying at a distance 29 from the plug, and generating either 

: +9 2836 i : 

an E,- or an H,-wave, R varies as sin? “yr where ’ is the wave length in the 
guide. These results are true, of course, for a guide of any cross section, and 
for any assumed current distribution. 





z 
i 
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4. Calculation of Reactance 
The e.m.f. method of calculating reactance uses the formula (3) 


(16) 24a ee # I(s)E(s)ds,* 

where the integral is extended along the antenna, and E(s) denotes the electric 
field along the surface of the antenna calculated for a current concentrated 
on the axis of the antenna. This procedure can be justified when the gap is 
along the length of the antenna, but it seems rather questionable to the 
writer whether it’ holds when the gap is across the antenna, as is usually the 
case in a wave guide (for the resistance, which can be calculated by means of 
the distant field, the procedure is less open to doubt). Assuming, however, 
that (16) holds, the reactance can evidently be calculated by means of the 
formulae in Section 2 giving the field. But as there will now be an infinite 
number of terms surviving in the series (7), the result will be expressed by 
means of an infinite series. Formula (16) assumes that the field is symmetrical 
about the axis of the antenna. This is evidently approximately the case if 
the radius of the antenna is small compared with its distance from the walls 
of the guide (except at the point of entry), and if it enters the guide normally 
(taking account of the image effect in the tangent plane). Otherwise an average 
of E(s) over the cross section of the antenna must be used in (16). 

We shall not give an explicit formula for reactance in the general case. 
Without detailed calculation, however, we can draw a rather interesting 
conclusion regarding the way the reactance of a transverse antenna, generating 
only an E,- or H,-wave, varies with the distance of the antenna from the plug, 
other factors being kept constant. Let the axis of the antenna be any curve 
lying in the plane z = zo, and let us take the integral in (16) along a parallel 
curve lying in the plane zs = zo + 7, where r is the radius of the antenna. 
Then use of (7), (11), (13) shows that the integral in (16) is, in its dependence 
on 20, 7, of the form 


(17) 2 [A meitm (220 +7) + Bret + Cet nl(20o+r) 4 Deu), 


where A, ,..., Dm are independent of 29, r. In the first and third terms 
of the summand in (17), r can be neglected. 

Suppose now, for definiteness, that only the E,-wave is generated. Then 
u, is real, but all the other u’s and ail the U’s (except possibly those for which 
the corresponding coefficient A,»,..., Dm is zero) are pure imaginary. 
Then (17) can be written 


(18) Angin + 2 Ane ttalm + 2 Cae-21Unie + D, 


where A,,..., D are independent of zo. The 2nd and 3rd terms in (18) 
will evidently be small compared with the other terms unless 2 is small, 


* Formula (16) is a little different from that given by Labus, but is more compact and is easily 
deducible from his formula. The negative sign appears here before iX on account of the time 
factor being taken as e- it, 
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i.e., unless the antenna is near the plug (when the above method would, in 
any case, break down). Excluding this case, therefore, we see from (16) 
that X is approximately of the form 


X = Ci cos (2un% + €) + C2, 


where C,, C2, € are independent of z9. An exactly similar result is obtained 
if only the H,,-wave is generated, with U, in place of u, . 

Hence if the distance from the plug to the antenna varies, other factors 
being kept constant, the reactance should vary approximately sinusoidally 
(except near the plug), the period being half the wave length in the guide. 
This is the same behavior as in the case of the resistance (see end of Section 3), 
though the phase will of course be different. Some experiments of Chipman 
and Burr (1) support this conclusion for the case of a straight transverse 
antenna in a rectangular guide. 
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PART III. THE RESISTANCES OF ANTENNAE OF VARIOUS SHAPES 
AND POSITIONS IN RECTANGULAR AND CIRCULAR WAVE GUIDES 


By W. Z. CHIEN 


1. Introduction and Summary 


The radiation resistance of an antenna inside a wave guide is defined as 
the ratio of the total power radiated to the mean square value of the feeding 
current. If the current distribution along the antenna is given, the first step 
is to calculate the field at the far end of the guide. When this has been done, 
the calculation of the total power radiated is just a straightforward integration. 

In the following section, we give a list of formulae for the resistance of 
various antennae inside rectangular and circular wave guides. These formulae 
are divided into four general cases according to the cross section of the wave 
guide and the type of wave transmitted. For each general case, we have a 
general formula for the resistance in terms of integrals depending on the 
shape and position of the antenna and its current distribution. The general 
formula (1), quoted from Part I, gives the resistance of an arbitrary antenna 
in a rectangular wave guide transmitting the Hi) wave only. The formulae 
(2), (3), (4) give the resistance of various antennae in circular wave 
guides. These formulae are calculated from the field distributions given by 
H. Buchholz (1). The purpose of keeping these general formulae in the list is 
that they offer a simple means of calculating those special cases missing in 
this list without knowing the actual field distribution. It should be noted 
that these formulae are but special forms of two more general formulae, which 
have been given by A. F. Stevenson in Part II. 


SST PET SE 





i 
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2. Results 


The Resistance of an Arbitrary Antenna in a Rectangular Wave Guide Trans- 
mitting only the Hi) Wave 


General formula: 
(1) R 


where 


(la) g = | , sin (= 20) cos (z 2) dy, (taken along the antenna). 
0 


2 
¥ 2407p é 


aly hms, 


The conditions for transmitting the Hijo wave only in a rectangular wave 
guide are in general: 
(1b) 


QS 


Hin 
Rio 


gi— 
A 

ol 

glam 


The notation is as follows: 
x, y, 2 = the rectangular co-ordinate system, origin at the center of the 
plug face, 
a the longer side of the guide along the x-axis. 
b = the shorter side of the guide along the y-axis, 
the wave length, 
y = 1 _ (x): “ y = the phase constant of the Ho wave. 
Xo, Yo, 39 = the co-ordinates of a point on the antenna, 
J = the current along the antenna, 
Io = current at the feeding point of the antenna. 
Special cases: 
la. A straight antenna of length 7 parallel to the y-axis, running from 
y = —b/2toy =1 — b/2 (Fig. 1). Assuming sinusoidal current 
distribution, 


_ 60d? E (z x» ) sin Ge Z ) tan (¥)] ohms 
~ maby ee tr ape 


where x» = the off-center distance from the y-axis to the antenna, 





Zo = the distance from the plug to the antenna. 


1b. The same antenna as in case la. Assuming uniform current distri- 
bution, 


240F | (z ) 7 i 
= —— | cos{— xo} sin |—,— 20) | ohms. 
aby a A 


1c. A circular loop, the axis of which is parallel to the x-axis, radius ro, 
center at (Xo, Yo, 20) (see Fig. 2). 


Assuming uniform current distribution, 


— 60ri— (Ty ary . (257 )} 
_— aby [cos (= 2) cos ( > 2) Ji Xr To ohms. 
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1d. A circular loop, the axis of which is parallel to the z-axis, radius ro, 
center at (xo, Yo, 20) (see Fig. 3). 


7 
| 
| 
| 
| 
| 
| 
| 





| 
| 
| 
| 
| 
| 


Z=Zy SECTION 


I, 
| 
| 
a a 
=O e-2T 
UNIFORM CURRENT DISTRIBUTION 
Fic. 3. 


Assuming uniform current distribution, 


( 3y2 2 
j= Rc sin (2 2) sin Ce Zo) Ji (z ro ] ohms. 
aby a a 


le. An American antenna of circular shape in the yz plane, being 
divided into four sections, radius ro and center at (xo, Yo, 20) (see 
Fig. 4). Assuming sinusoidal current distribution in each section so 
that mro = X, 


_ 61440mr — (7. 2 P 
R= “a [cos (2 2) cos ( - 30) Fe) | ohms, 





where 


a 
dx? 
if. An American antenna of the circular shape in the XY plane, being 

divided into four sections, radius ro and center at (xo, Yo, 20) (see 


Fig. 5). Assuming sinusoidal current distribution in each section so 
that mro = X, 


_ 614407, . (ar ~~ A\?? 
R= ~~ [ sin (2 xo) sin ( x v0) F ()] ohms. 


F(x) 





» Wo = Xr. 














g 
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The resistance of an Antenna Parallel to the Axis in a Circular Wave Guide 
Transmitting the Eo, Wave Only 


General Formula: 


30 ? 
@) R= sR Haat LPo7Ln Gt) | ohms 


where 
3 3 
(2a) the numerical factor = FRG = 11.28. 


(2d) go = | 7 cos (= Yate) dzo (taken along the antenna). 
0 


The conditions for transmitting the Eo: wave only in a circular wave guide are: 


(2c) antenna parallel to the axis, 


2 
(2d) ju = 2.405 <: > < 5.520 = jos. 
The notation is as follows: 
a = radius of the wave guide, 


ju = 2.405 = the first root of Jo(x) , 


= |i — (GAY = oft — (0.38272Y 
a1 = 4/1 ) = 1 (0.3827 *) ’ 


27 4 

— 1: = the phase constant of the Eo: wave, 

ro, 29 = the co-ordinates of the antenna, ro = the off-center distance, 
Zo = the distance from the plug, 
I = current along the antenna, 
Io = current at the feeding point. 


Special cases: 
2a. A straight antenna of length / parallel to the z-axis, running from z = 0 
toz =/. Assuming sinusoidal current distribution (Fig. 6), 


308.0 sin K (1 + =o) sin iF a —- re) | oi 
R = —_ Sy Jo . ohms. 


™ sin (22) 


480 
: —, = 308.0. 
li Ji(jm)}? | 


2b. The same antenna as in the case 2a. Assuming uniform current distri- 


bution, 
At 2 ‘ 
R = 1.652 - a iyi ; sin (= ral) Jo (ie a) ohms. 


15 Jor - | 
—~| = 1.6852. 
ie & Ji(jor) 5] : 


9 
” 











Ct RY 
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DISTRIBUTION DISTRIBUTION 
Fic. 6. 


The Resistance of a Circular Loop Antenna Coaxially Placed in a Circular 
Wave Guide Transmitting the Hy, Wave Only, Assuming Uniform Current 


Distribution (Fig. 7) 


4807? 6 ro 


eo ee ee , 
©) R= Faciay ent L¥ (Sa) sin (FF oe) J ot, 





I, 
| I, 
6:0 @:27 6:0 @:2r 
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where the numerical factor is 

4807? 
J6 (Jor) 
The conditions for transmitting the Ho wave only in a circular waveguide are: 
(3b) symmetrically placed loop antenna with uniform current distribution, 


(3a) = 2.921 x 104 


(3c) ju = 3.832 < ore < 7.016 = jis. 


The notation is as follows: 
jo = 3.832 = the first root of Jo(x), 


ro = the radius of the loop, 
a = the radius of the wave guide, 
maa 5 8, paren re aes 
ono = aii — (2%) =,/1 - (0.6098 2 ‘ 
You V G des) Vv a 
2 F 
> Yo. = the phase constant of the Ho wave. 


The Resistance of an Arbitrary Antenna in a Circular Wave Guide Transmitting 
the Hy, Wave Only 


General formula 


Sd cee i ) : 
“7 R = GES D ya (er ~ Ge) ohms 
where 
240 
4a) the numerical factor = —>.——+5———— = 296.6, 
(4a i sa Ji(ju) (ait — 1) 


= E ae . (24, \ dro (taken along the 
o t= / To” (2) cos Oy sim (= Yuto) 7 antenns), 


sl cg hae ate » fle 
(4c) Yip = [i J; (ii: z2) sin 4 sin (= viute)ro dO). 


The conditions for transmitting the 17;; wave only in a circular wave guide are: 


(4d) for an antenna having a component in the axial direction: 


Ss 27a ; 
Ju = 1.841 “> < 2.405 = Jes 

(4e) for an antenna having no component in the axial direction: 
fe 27a ea a 
jin = 1.841 < — < 3.054 = 793 2 

The notation is as follows: 

a = the radius of the wave guide, 
ro, 00, 20 = the co-ordinates of a point on the antenna, 
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ju = 1.841 = the first root of J; (x), 


Fad ~ {i a = S 
Vu 4! (25) RP (0.2930 *) ; 
ai Yi1 = the phase constant of the Hi, wave. 


Special cases: 


4a. A straight antenna of length / placed along a diameter of the wave guide 
with one end at the surface of the guide. Assuming uniform current 
distribution (Fig. 8), 








UNIFORM CURRENT 
DISTRIBUTION 


Fic. 8. 


et fe ? s in\\? 
R = 296.6 v, [ sin y Yn | { Bin) —B (( - ni)| ohms, 
where 


B(x) = Jo(x) + | Jo(a) da, B(ji,;) = 1.382. 
0 
, 1 athe cs ; 
For the special case ] = 2a: R = 2259 yi, [ sin (< viuto) | ohms. 
il 


4b. A circular loop of radius ro placed coaxially in the guide at a distance 20 
from the plug. Assuming sinusoidal current distribution (Fig. 7), 


4 
Re 14508 = 
Ory ii (1 47° rn) 
— 0 





r2 


[ sin( 5 als (= A er 
sin r 71120 an r:” 1 ING onms, 
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where the numerical factor is 


2 
240 jij 167? = 1.588 X 105. 


(Ai — 1)J7 (ju) 
Reference 
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PART IV. THE IMPEDANCE OF A RECTANGULAR WAVE 
GUIDE WITH A THIN ANTENNA 


By L. INFELD 


WitTH APPENDICES BY J. R. POUNDER AND A. F. STEVENSON 


1. Introduction 


A wave guide of rectangular section extends to infinity in one direction and 
is closed in the other direction by a plug. The guide and the plug are assumed 
to be perfect conductors. 

A thin, cylindrical antenna is placed inside the guide. To make our argu- 
ment simple, we place this antenna parallel to the plug and parallel to the two 
other (shorter) walls of the guide, midway between them. Fig. 1 shows the 
choice of our co-ordinate system. In it, the lower end of the antenna has the 
co-ordinates (0, 0, 0) and the upper end has the co-ordinates (0, J, 0) where / 
is the length of the antenna. The z-axis is perpendicular to the plug, and 
points toward the open end. 


| 
| 
| 
| 
=) 


—— —— m— — —— - 






x 
Fic. 1. 
The parameters characterizing our problem are: 
a,b, the dimensions of the guide, 
2a 
k= - (\ = wave length), 
1 = the length of the antenna, 
Po = the radius of the antenna, (we assume pok < 1), 
Zo = the distance between the plug and the antenna. 
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Our problem is to find the impedance of this wave guide under the assump- 
tion that the current distribution is sinusoidal. In this respect this paper is 
a continuation of Part J, in which the problem of the resistance of a rectangular 
wave guide was solved. However, the method of finding the impedance, and 
therefore both the resistance and the reactance, has turned out to be so 
simple that I intend to formulate this part independently, briefly repeating 
the reasoning concerning the images. 

Our particular assumption about the position of the antenna is not essential. 
It would be easy to generalize our argument in this respect, as was done in 
Part I for the case of the resistance (see Appendix A). However, the assump- 
tion that the antenna is of cylindrical shape and that the current distribution 
is given is essential for our argument. Again, the assumption that the current 
distribution is sinusoidal is not necessary and the argument in the case of a 
constant current distribution and small antenna would be even simpler (see 
Appendix B). 

Some of the results found here are valid for a wave guide of arbitrary 
shape as long as the antenna is parallel to the plug. This is shown by Steven- 
son in Appendix C. Especially the conclusions formulated in Section 8 are 
general and may be of interest to someone who does not care about the detailed 
calculations. Therefore, I tried to write Section 8 so that its contents may be 
understood without reading the previous parts. 


2. The Lattice of Dipoles 

We use the method of images. The image of a dipole with respect to a 
perfectly conducting plane is obtained by taking the geometrical image and 
then reversing the direction. The dipole and its image give a field satisfying 
the condition 

sae 
(2.1) E, = 0 
+ 
over the plane, £, being the tangential component of the electric field. 

Let us first consider the case where the rectangular wave guide has no plug, 
so that it extends to infinity in both directions. Suppose there is a dipole, 
inside the guide, placed along the y-axis. It is easy to find a doubly infinite 
system of images (all external to the guide) such that the lattice consisting 
of the dipole and its images gives a field satisfying the condition (2.1) over the 
walls of the guide. If we now add a plug, the condition (2.1) will be satisfied 
over it if we add to the above lattice its image in the plane of the plug. 

We place a dipole with components (0, A,, 0) at (0, yo, 0). Fig. 2 shows 
part of the lattice in the plane z = 0. 

If we let m and n take all integral values, positive, zero, and negative, we 
may represent the components and positions of all the dipoles in the plane 
z = 0 as follows: 

(04,0) at (2ma, +y¥o + 2nb,0) —_ two basic lattices, 
(2.2) 
(0, —A,,0) at (2m +14, +y¥0 + 2nb, 0)—two basic lattices. 








2 
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The images with respect to the plane of the plug are given by this table 
by changing the sign of A, and replacing z = 0 by zs = —2z9. Thus we have 
eight “‘basic’’ lattices determining the field in the guide. 





Fic. 2. 


3. The Impedance of a Free Antenna 

Before we return to our subject of the wave guide, we shall derive an expres- 
sion for the impedance of a free antenna driven at the lower end. We assume 
the antenna to be as in Fig. 1, but we shall ignore the wave guide. 

We shall now use Sommerfeld’s formula (4), which allows us to write the 
Hertz vector (in Heaviside units) of a dipole at (0, yo, 0)in the form of an 
integral :* 

a +o 
@@ tf . aos 
u vy ,, 9 > 
ee ele ad \ eiay g-iavo FI Jk? — of da, 
a Sar a (pv ) 


-« 


(3.1) 
P=ae+2, 7 = p+ (y — yo)’, 
VE — oe = positive imaginary for a? > ’, 
Hy = Hankel’s function of the first kind and order zero. 
We omit here and later the time-factor e~**. 
If we have, instead of a dipole, a finite antenna, then the Hertz vector of 
this antenna can be found by replacing A, in (3.1) by 


i sin kd — yo) 4, 
ke’* sink) °>°" 

and integrating over dyo from 0 tol. Here Jy is the current at the driving 
point (yo = 0), and the expression (3.2) means the assumption of a sinusoidal 


current distribution. 


(3.2) 


* Our formula follows from Equation (3) on p. 416 in Watson's book by putting wh = —}4, 
v= 0, 
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It is convenient to introduce 
I 
(3.3) f(a) = [eo sin k(l — ¢)d¢ 


0 


' 
: 
& 


ioe <alnepinerint seca NI oc i acyl 


-@) | 


Using (3.3), we have for the Hertz vector ®, of the finite antenna: 


os So jay FT eee 
(3.4) ig gee |r ee HY (pk? — a2) da. 
The next step is to introduce the y-component of the electric field, the only 
one in which we are interested: 
. ‘ ah 
Tw = - ®, + “ay? 
(3.5) 


—I : . i aaa 
= 8rck a fw = a’) f(a) oC HY (pV k? = a’) da. 


Our last step is to calculate the impedance (Zoo) of this free antenna by 
the e.m.f. method: 


l 1 
Y = x. , ae ae i, i om 
(3.6) Zo = — ia Bao Idy = aa [20 sin k(l — y)dy. 
From (3.5), (3.6), (3.3), we obtain: 
9 pela ae ci " Pi bLS ae (1) 1a ae 
(3.7) Zoo = Sack seam | f(@) (k a’) Hy (poVk a’) da, 


where f(q@) is the complex conjugate of f(a). 


Equation (3.7) represents the impedance of a free antenna. 


4. The Mutual Impedance of Two Antennae 


Before we return to the wave guide we must solve the problem of the 
mutual impedance of two antennae. 


Let us assume that the two antennae are parallel to each other and of the 
same length 7. Let p be the distance between them. Then: 


l 
: 1 
(4.1) Zap = — ces (p) I dy. 
0 


Here Z,4z is the contribution to the impedance of antenna (A) coming from 
antenna (B). The indices (A) and (B) denote the current and the field 


i 





; 


eR 
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belonging to antennae (A) and (B) respectively. Finally we may remark 
that our formula assumes po < p and therefore it is immaterial whether p 
means the distance between the axes of these antennae or between lines 
parallel to them and lying on their surfaces. 

We shall now consider two cases: 

Case a: Each point of antenna (B) is obtained from a point of antenna (A) 
by replacing 
(4.2) 0,¥0,0 by do, ¥o + bo, 60, 

ao, bo, co, being constants. 

The calculation of Z{%% is very similar to that of Zo. Where we had 
y — yo before, we shall now have y — yo — bo. Secondly, where we had po 
before, we shall now have 
(4.3) p= Vari t+ce. 

Thus, instead of (3.7), we now have: 

1 bigs iia ceca 
4.4) Z%, = ————— | ei” f(a) f(a) (#2 — a?) HS? (pV — a2) da. 
(4.4) AB Suck oint ul F(a) f(a) ( ) Hy” (pv ) 

Case b: Each point of antenna (B) is obtained from a point of autenna (A) 

by replacing 
0,970,090 by ao, —yo + bo, ce. 


The calculation of Z{, will differ from that of Z{% in one respect: in the 


integration with respect to yo we shall have 
I 
(4.5) ei sin R(l — yo) dyo = f(a), 


0 
instead of 


l 
[oo sin kil _ Yo) dyo = f(a) 
0 


as before. 


Thus we have for Z°}, : 


Zo) 1 os f(r) )2 (2 2 ( a 
(4.6) ies = aan | (f(a) )* (2 — a2) HY (pV — a2) da. 


@ 


Both results, (4.4) and (4.6), will be needed in our later argument. 


5. The Impedance of a Wave Guide 
We summarize the result of our previous section: The mutual impedance of 
two “‘parallel’’, identical antennae (case a) is: 
+a 
(5.1) wa = | e'% o(p,a) da. 


—o 
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The mutual impedance of two ‘‘antiparallel’’, identical antennae (case 0) is: 


(5.2) iB - [em V(p, a) da. 


eo 


Here we have, according to (4.4) and (4.6): 


g(p, @) = ga fla) fla) (kt — a) He (pV — a), 
(5.3) 
Vp, a) = ge (Fla))? (# — a) HY (pV = 0), 


and according to (3.3) 


: ei — cos kl + i@ sin kl 


(5.4) f(a) = ~—7 
~~ (GY) 1 


The self-impedance of an antenna is obtained from (5.1) by putting there: 
(5.5) bo = 0, p = Po. 

Now we know that the impedance of a wave guide is equivalent to the 
impedance of an antenna surrounded by eight basic lattices. In (2.2) we 
wrote down the co-ordinates of the dipoles of four basic lattices. If the 
dipoles become finite antennae, the four antenna-lattices consist of ‘‘parallel”’ 
or “antiparallel” antennae. The remaining four lattices are opposite in 
direction and placed in the plane z = —2z9. Thus we have, because of (5.1), 
(5.3), and (2.2) for the impedance of an antenna in a wave guide: 


m=-0 n=-© 


lin Df fo [ 2(p0 + 2ma, a) + W(po + 2ma, a) 


(5.6) — ¢( (2ma + 1)a, a) — Y( (2m + 1)a, a) 
— eV Omar + Baal a) — WV may Cay, «) 
+ eV (m+ 1)aP+ 2z0)% @) + WW (m+ 1)a) + (220), aw) | dew 


We have here, in square brackets, eight expressions, each corresponding to 
one of the basic lattices. The terms with po, for m = 0, give the self-reactance, 
and in them pp must not be omitted. Otherwise, for m # 0, it is immaterial 
whether po is omitted or not. 

Our expression for the impedance, (5.6), though representing the solution 
of our problem, is too complicated to be used in actual calculations. Our aim 
is to change (5.6) in such a way that it will be applicable for practical purposes. 
We shall do it in detail for the first expression: 


+a 
(5.7) >) [= ¢(po + 2ma, a) da. 


@ 





EF DETER RE ae 
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We introduce a new variable 


b 
go =-—a. 
w 

Then instead of (5.7) we have 
(5.8) 7% % ¢ (6 + 2ma, = c) do. 

We can now apply Poisson’s formula (1), replacing (5.8) by 
(5.9) FS 3, (m0 + 2ma, ). 

For m ¥ 0, we have p = V/4m?a? = 2|m|a. Therefore, the summation 
from m = — © to m = + © can be replaced by a summation from m = 0 
tom = ~. If we look back at (5.3), we see that g(—a@) = g(a); therefore 
g(— nn) = g(n). Again the summation from n = —o@ to m = @ can be 


replaced by summation from 0 to ». Thus (5.9) takes the form 


> bm On ( po + 2ma, w) , 


On 


Tv @o 
(5.10) 3 > 


where 


0 
6, = Zior's x 0, Gy = I. 


We have to repeat the same reasoning for the remaining seven expressions 
in (5.6), each time using Poisson’s formula. Having done this and having 
introduced for g and w expressions (5.3) and (5.4), we obtain for the total 


impedance: auis 
1 oo (cos kl — cos =m) 
Z ee eS a em 7 Pe gig oe 
Ot) 4” gay ant oa8 vo Tae ne (3) - 
\bk 


HY c + ma) ve- (| 
— HS (ome (224)*) ( .e ))| 


This is the final result. If expressed in ohms, instead of Heaviside units, it 
must be multiplied by 1207. 





6. Discussion of (5.11) 
The great virtue of (5.11) is its quick convergence. Let us assume 


TH 





- (6.1) =a” y> 1, mak = p>>1. 
Then 
Ee 
6.2 HY (ipv)| ~~ 
(6.2) | b(t )| V/ wy 


which shows that the evaluation of (5.11) will be comparatively easy tohandle. 
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Let us imagine a point lattice (Fig. 3), where each point corresponds to one 


expression in the sum (5.11). 


IMAGINARY CONTRIBUTIONS 
e 
° 


the wave guide are 


Gi 


\' - (®) - 


(6.3) 


vi=( 


Let us further assume that the dimensions of 


. . e e. . 
NEGLECTEO 


CONTRIBUTIONS 


OMPLEX CONTRIBUTIONS 


De 





Fic. 3. 


such that it lets through only the Ho, wave, and that 


Am\? is imaginary for m = 2,3,----, 
2a and real for m = 1, 
An\? . . : 

— (55) 38 imaginary fOr Gee 1D. 259 


Because /7§” of an imaginary argument is imaginary, we see that the con- 


tribution to the impedance is purely imaginary if nm 2 1. 
along the m-axis give a complex contribution. 


Only the points 
Because of (6.2), the most 


important contribution will come from the points on the m- and n-axes. It 


was found that all others can be ignored in practical applications. 


Therefore, 


we shall now focus our attention upon the contributions to (5.11) from these 


points, which are as follows: 


(a) m = 0, n = 0: 


tan? — 


(6.4) 
“Abck 


- (I$? (pok) — H§? (2z0k)), 


e 
i 
' 
& 
f 
; 
ft 
' 





RP 





i 
& 
: 
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(6) m=0, n= 1,2,----: 


1 a (cos kl — cos amy as 
Aisin tee aie 
| (oot) (F 


2bck sin? kl z eo (Fz) 


bk ae 
(6.4) (Cont'd) = H19°( sok] — (Fz) | 
(c) me = 0; wm = 42.---=: 
tan? kl/2 2 ; 
a x (— 1)" [2§? (mak) — H§? (RV (ma)? + (220)?)] . 


As Hj” (x) behaves like log x for small x, the (a) part gives the logarithmic 
term which appears in all calculations of impedance. In the (0) part, we shall 
have to sum a few expressions only, if pok is of the order of, say, 1/10. Finally 
(c), that is the summation along the m-axis, is very simple, although it does 

; Ce cs 
not converge exponentially. However, even (c) converges as } 2’ Since 
for great m we can write . 

(— 1)"(H§? (mak) — HS? (kV (ma)? + (220)2)) ~ ein oa H (mak) 
ma sie 
a 
We shall see later how to deal with part (c) in actual calculations. 

If we are looking for reactance only, we shall have to replace H§? by i Yo in 
(6.4) (a) and (c). 

7. The Resistance 


The resistance is the real part of (6.4) (a) and (c) for po > 0. 
tan? ~ «wo 
(7.1) R= ce ae (Jo(mak) - Jo(k/ (ma)? + (220)*) Jer, 
where R denotes the resistance. 
We shall again use Poisson’s formula, this time replacing summation by 


integration. Then we have: 


» Rl 
tan* — ? 


capa 
4bck h=-© 


48 


(7.2) R= e-ita(h- veal Jo(aak) 


— Jo(k-V (aa) )® + (250)*) de 

We can replace the Jo's by H§”’s and calculate R as the real part of 
» kl 

an = « : 


a pe enitm(h-1 al 715? (ak) 


(7.3) 


— HY (kV (aay + Be) Jace. 
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Now we can use Sommerfeld’s formula (3.1) with 
: 2120 





yo — 2r(h — 1/2), p — tak, y > 0, R > 


We have, according to this formula: 


eo 


Sone a rT IE eee mer me peeR 


+o +a 
[ersreveminy eats —_ [eve 1a bee (kv/ (aa)? + (220)2) dex 


@ 


(7.5) 


-230  ——$—_—_—$—__—_——— 
2 2 q V~ ak) +4012)? 
e 


~ ix/— (ak)? + 4? (h — 1/2)? i\/— (ak)? + 40h — 1/2)? 
We see that this expression will be real only if 
(7.6) =O ‘or “hk => 1. 
In this case we have for the above integrals (for 4 = 0 and h = 1): 


: & Mi 2: 
in gg __ #008 (22k) - (5) 








The result is 


2 tan? - sin? (cst, ( = )) 


( ve 8) R = SS ee 
9 i 
abck 1 (=) 


which is in perfect agreement with the result obtained in Part I. 


a 
=| a 





8. The Impedance Circle 


We shall present our reasoning here in the form of theorems, so that the 
reader interested in their formulation only can omit the proofs. 


We assume the dimensions of the wave guide (a, 6) and the wave length (A) 
and the thickness of the antenna (2p 9) as known and kept constant. The 
characteristic impedance of the line (Zo) is also known. Our problem is to 


find J) (the matching length of the antenna) and 2 (the distance of the antenna \ 
from the plug) so that the resistance will be that of the line and the reactance ( 
will be zero. It is convenient to take as the unit of impedance that of the 

line (Zo). In this unit our problem is to find Jy , 29 , such that 

(8.1) Rl, X= O,whereZ = R — 4X: ( 


This we shall refer to as the matching problem. 


Theorem I 


The resistance can be expressed in the form 


(8.2) R = Ail) - A(l) cos (z v9); Ww 
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where 


ect (%) 
‘isi | A), = ey tan x)” 
( 


The proof follows immediately from (7.8), if we assume that the 
impedence of the line (Z») is expressed in ohms. 
Theorem II 


If zo is not small compared with \, then we can express the reactance in the 


form 
(8.4) X = A’(l) + A(l) sin ( 2) , toa very good approximation, 


where A’(/) is a certain complicated function given explicitly in the course of 
the proof. But—this is essential—neither A nor A’ depends on Zo . 


Proof 

The reactance X is the imaginary part of —Z. From (5.11) we see that 
—Z can be split into two parts, one which does not depend on 2 and the 
other which does. The part of Z which does not depend on 2» is 


(cos = — cos al) 


(— 1)" bn 5, - a (s ~~ xX 
FE) 


AY io + ma), |i! = (7) | ; 


It was shown in Section 7 that the real part of this is A, as defined in (8.3). 


(8.5) =." © 


Zobk sin? Rl 2 t20 





We shall therefore denote the expression (8.5) by 


A —iA’, 
where A’ is real. Thus we can write: 
(8.6) X = A’ — Im(C), 
307 a (cos ee =) 
(87) = Fob sin? Bl ey 6 1" Om On ae a 
bk 


Hs» [ env (ma)? + (20)? ] ; 


: - rn\? 
where k2 = BR - (=") : 
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Before, in Section 7, we calculated the real part of C. We shall now cal- 
culate in a very similar manner the imaginary part of C. 


Using Poisson’s formula (1) we can write for C 


(8.8) Pa (ws kl — cos a) “ = 
, 307 x4 anneecioicinieningeenien: ae en i2m(h-1/2)a Xx 
~~ Zane sin? ki ,o oe ~ hom 
bk il 


I Hi [ka (ma)? + (220)? da. 


Now we proceed exactly as in Section 7, (7.4), (7.5), using Sommerfeld’s 
formula (4). Therefore: 


(8.9) sa 7 (cos re om) - 4 ee 
i in? setae re MO Ie pagan se ita ea ide 
bkZo sin® kl =o ms (3) hao in/ — (akn)® + 40*(h—1/2)? 
bk 


For n > 0 the exponent of e is real and negative. If the antenna is not too 
near the plug, then the only contributions to C which have to be taken into 
account are those coming from = 0: 
60m tan? - 2 Vani 
eee: 5 ee ee ee ee 

bkZo — nz—0 in/ — (ak)® + 40°%(h — 1/2)? 


For the same reason we can neglect all expressions in this sum with the 
exception of those for which 


(8.10) C= 











(8.11) kh = 0,.0r & = 1. 
Therefore: 

(8.12) Im(C) = — Asin (277), 

and finally 

(8.13) X=A'+A sin (* me) 


which proves our theorem. 


Theorem III 


If we change 2 and keep all other parameters (a, b, A, po, 1) constant, 
then the points characterized by the corresponding pairs (R, X) lie on a circle, 
called the impedance circle. This circle is tangent to the X (reactance)-axis 
and the co-ordinates of its center are (A, A’). 


This can easily be seen from (8.2) and (8.4). Indeed, eliminating 2) from 
these equations, we have 


(8.14) (R — A)? + (X — A’)? = A’, 


which proves our theorem. (See Fig. 4). 
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What is the physical significance of the symbols appearing in (8.14)? Rand 
X are the resistance and reactance corresponding to given values of / and Zo. 
As 20 changes, the point (R, X) moves along the impedance circle. 2A is the 


x 


Fic. 4. 
maximum resistance for a given / and variable 25. A’ can be interpreted as 
the reactance at the point for which the resistance vanishes; indeed for 


2ry 
A 


we have R = 0,X = A’. 





3 = se; s = 1,2,3,----, 


9. The Matching Problem 


We can now turn toward the solution of the matching problem. We have 
to find a circle, and a point on the circle, for which : 


(9.1) R=i, 4 = € 
Putting these values into (8.14) we have 

Am 
C = =, on 6 
(9.2) A 2 + 3 


Therefore the center of the impedance circle must lie on the parabola (9.2) 
for the matching problem to be solved (Fig. 5). The remaining problem is: 
how to find the right point on the parabola? Theoretically speaking, we 
ought to regard (9.2) as an equation for the determination of J) , as A and A’ 
are known functions of 1. 

If we look at (9.2) we see that the difficulty is to find A’ as a function of /, 
as the other expression, A, is very simple. But A’ isa part of Z; more explicitly 
the imaginary part of that part of —Z which does not depend on 2. Let 
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us rewrite from (6.4) the part of —Z in which we are interested, that is, the 
part which does not depend on zo. This is: 








a 
ae 
| (a) ———— Yo( pok) 
‘i 
93) | | cor @ (cos A — cos :¥) — 
| (>) Fez, sin? bt ve A) 118°(0 oka] ~— (F) ) 
| + 





| 607 tan? - x 
| (c) gage a (— 1)™ Yo(mak). 








Fic. 5. 


Thus we have here the explicit expression for A’ which was used in all 
further calculations. In computing, the only difficulty encountered was to 
calculate (c): 


(9.4) x (— 1)™ Yo(mak), 
m=1 


which converges very slowly. This mathematical problem is treated else- 
where (2), where a table and graph of (9.4) is given. 
To make our procedure clear we shall take the following special example: 
(9.5) a 7.213 cm., 6 = 3.404cm., A = 10.0 cm,, 
po = 0.397 cm., Zo = 80 ohms. 
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Then A’ is a function of / only. Putting different values of 1, we have 
4(1 + A’(l)?) as represented by Fig. 6. Next we draw A (which is much 


0.56 


0.55 


0.54 





0.53 












0.52 





2.00 1— 25 2.50 


——scm. 


Fic. 6. A and 4(1 + A”) as functions of l.X = 10.0 cm., a = 7.213 cm., b = 3 404 cm., 
p = 0.397 cm, 
simpler) as a function of /. The intersection of these two curves has the 
co-ordinate J) , satisfying 


(9.6) $(1 + A’(lo)*) = A(lo). 


Once we have found /) , we can easily find the set of go's from (8.2) and (8.3): 


30 le wh _ 4ryZy 
R=1= thevd, tan x (1 cos X ) 
Thus the method allows us to find /) , 9, for which R = 1,X = 0. In our 


particular case (9.5) we have 


lo = 2.316; 2 = 9.839 + As . 


- gs = 0,1,2,---.- 
27 


In this way Jo, 2 were calculated for two wave guides, for different X's, 
po's, Zo's. The calculations were made by Chien and Pounder. It proved 
that the omission of all the terms for which m, n > 0 was justified. In the 
expression (b) of (9.3) it was sufficient to take some six terms. 
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APPENDIX A 


Generalization of the Impedance Formula for an Unsymmetrical 
Position of the Antenna Relative to the Sides of the Guide 


By J. R. POUNDER 
We assume now that the antenna is placed at a distance xo from the center 


of the guide (Fig. 7). Let us generalize the previous reasoning step by step. 
We have four lattices of dipoles in the plane z = 0, 


0, A,, 0 at xo + 2ma, +¥0 + 2nb, 0 





(A1) 
0, —A,, 0 at 2m + 1a —%X0, + Vo + 2nd, 0 
Y 
F 
| 
| 
4 
| 
| 
| 
| 
x 
Fic. 7. 
and four other lattices in the plane z = —2z9 obtained by reflection in the 


plug and by changing A, into —A,. (Fig. 8) 
The impedance of the antenna, corresponding to formula (5.6), is now 
by exactly the same method as before: 


Z= > >> / e—i2nba dal o( po + 2ma, a) + (po + 2ma, a) 


— ¢g(po + 2m+1a — 2x0, a) — W(po + 2m + 1a — 2x, a) 
(A2) = @(-V(2ma)* + (20)*, a) — WCW (2ma)® + (220)* , a) 

+ g(Ni2m +1 a — 2x0)? + (220)? , a) 

+W(N@Q@m + 1 a — 2x0)? + (2e0)? a) | - 











: 
i 
F 
f 
I 
; 
' 
' 
; 
i 
; 


en 


——— 


EE 


SERRE EEE 





CHIEN ET AL.: CONTRIBUTIONS TO THE THEORY OF WAVE GUIDES 125 





Fic. 8. 


Again by the use of Poisson's formula (1) we obtain: 


30m es ra (cos kl — cos my 
bk sin? kl pt note ' ae ay a 


bk 
(rem = 
_ HY (os + Im Fa — 2x) le — (™ (7 y] 
— HY [ce (2ma)* + (224) )(# nm F # - (7))] 
ine [ a (Gm FA a — 2)* + (250) ) (# - ())]} : 


As can be shown by computation very similar to that of Section 7, but 
somewhat more complicated, the real part of (A3) is: 


al 9 [ SSe\ .... fanyae\ _ — 2720) 
(A5) R = 2A cos (2 ) sin (A ) = 2B sin ( x , 


where A is defined by (8.3), and 





Z (inohms) = 


(A3) 








TXo 


(A6) B = A cos? 


This formula agrees with the corresponding one in Part I. It is obvious 


_ that for x9 = 0 this expression for the resistance reduces to that given in (7.8). 


It can be shown (with the same assumption as before, namely, that zo is not 
small compared with \) that the reactance can be written: 


(A7) X = B’ + Bsin ary where 
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He? (a + 2ma)k V3 es (zy | 
— He | (o + 2m + 1a—2x)k J! _ (#2) |}. 


Now (AS) can be written 





ATryzo 30? TXo kl 
AC in — i. cei hail Tee ge oe 
(A9) R B (1 cos > ) ; hei cos tan 9 


It is evident from (A9) and (A7) that all that was said before about the 
impedance circle is true also in this more general case. 


APPENDIX B 


Antenna of Length / = b 





By J. R. PounDER 


In the case in which the length of the antenna = 6, one usually assumes 


uniform current distribution (see Fig. 9). The problem is not new, having 
been solved by Schelkunoff for a infinite guide without a plug (3). One of 
the two methods used by Schelkunoff seems to correspond to that presented in 


this paper. 





Fic. 9, 


The same argument may be used here as before, or a slightly different 
treatment may be given which replaces the separate antenna-images by 
The essential point of either method is that the 


infinitely long antennae. 
Here we shall outline the 


double-sum is replaced by a simple summation. 


argument of the first method. 
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Because of the assumed uniform current distribution, we obtain instead of 
formula (3.3), 
1 a e—iad 


(B1) fla) = —> 


1a 


By the same argument as before we arrive at formula (5.6), but the next 
step (the use of Poisson’s formula) gives us an infinite series corresponding 
to (5.11) in which all terms involving »#¥ © vanish. What remains is: 


Z == > { H3” [ (po + 2ma)k | 
(B2) — HS? [(2m + 14+ po — 2xo)k] — Hi? [kV (2ma)? + (220)? ] 





+ HY [k Vm +1 a — 2x0)? + (220)?]} . 


From this we can compute, just as before (using Poisson’s and Sommerfeld’s 
formulae), the resistance and the reactance. We obtain 


oe Ocoee FHL _ re SEY ae ee 
(B3) RK = a cos? — (1 cos > i y = 1 ia)? 


Es b 9 WXy . 4ryx% be 
(B4) X = ao ’ sin — - tee ea { H§” [ k(po + 2ma) | 


— Hg [k(po + 2m + 1a — 2x) ]}. 


If we consider only the special case in which x9 = 0, and at the same time 
change from Heaviside units to ohms, these expressions become: 


5) 2 me 120d fy _ ong 21%) . — Dos 
(B5) R= “ (1 cos > ) = Zo (v D cos x ). 


120mb .. 4myzo 60m SQ any 
~ sin —\ 7 ae (— 1)"Yo(Rk(po + ma) ) 


z (D' + Dsin its) , 


where Zo is the impedance of the line, and D 


(B6) X = 


Il 


12075 
ayZo 


Obviously, the theorems about the impedance circle hold as before. 
The matching problem can be stated as before: How shall we make 
(B7) R= Zo and AX = 0? 
We eliminate 2 as in Section 9, getting 
(B8) D’? +1 = 2D. 


We cannot now satisfy (B8) by choosing /, since / has been chosen as equal 
to b, but perhaps we can find some other parameter instead. Assuming a 
given wave guide, that is, given a, 6, and \, we may ask: 

Having adjusted 2) so that R = Zy), can we adjust po, the thickness of 
the antenna, so that X = 0, that is, so that (B8) is satisfied? 
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It turns out that it is quite easy to match the line by choosing po. Here 

we shall give the matching conditions for the guide described in (9.5), that is, 
A. = 10.0 cm., a = 7.213 cm., b = 3.404 cm. 


Po and 2, giving the matching conditions, are represented as functions of 
Z, by Fig. 10. As an example, imagine that po = 0.1985 cm.; then we see 
from Fig. 10 that the matching requires Z) = 60 ohms and 2) = 6.15 cm. 


0.25 6.5 


0.24 


0.23 











60 
Z)- OHMS — 


Fic. 10. po and zo as functions of Zo. X = 10.0 cm.,a = 7.213 cm., b = 3.404 cm. 


If experiment were to verify these results, it would mean that the current 
distribution is uniform. 


APPENDIX C 


Resistance and Reactance as Functions of Distance between Antenna 
and Plug for a Wave Guide of Arbitrary Cross Section 


By A. F. STEVENSON 


The results (8.2) and (8.4) for resistance and reactance can be written 
(C1) R = A(1 — cos 2uy) , 
(C2) X = A’+A sin 2uz, 


where Zp is the distance of the antenna from the plug, u = 27/)’, N’ being the 
wave length in the guide, and A, A’ are independent of 29. We wish to point 





i 
E 
f 
| 
j 
: 
| 
: 
; 
i 





CHIEN ET AL.: CONTRIBUTIONS TO THE THEORY OF WAVE GUIDES 129 


out that the results (C1), (C2) hold for a wave guide of arbitrary cross section 
and an antenna of any shape which lies wholly in a plane distant 29 from the 
plug, provided that the antenna is not too near the plug, and provided that 
the dimensions of the guide and manner of excitation are such that only a 
single E- or H-wave is transmitted. The current distribution in the antenna 
(provided it is independent of Zo) is immaterial. 


In Part II, it was, in fact, shown that, under the conditions stated, R is 
proportional to sin? uz, which is equivalent to (C1), but (C2) was only 
established in the less precise form 


X = Ci cos (2uz + €) + Ce. 


The result (C2) follows easily, however, from (16) and (18) of Part II. 
For, disregarding unimportant terms, these give for the impedance in the case 
considered 


(C3) R — iX = Berum+C, 
where B, C are (real or complex) constants independent of 29. Put 
(C4) B = B, + 7iB., C=z=C+ Ce. 


Then equating the real parts of each side of (C3) we obtain 

R = B, cos 2uz) — Bz sin 2uzo + Ci. 
Comparing this with (C1) (already established), we see that we must have 
(CS) C. = -—B, =A, Bz =0. 


Substituting (C5), (C4) in (C3), and equating the imaginary parts of each 
side, we obtain 

A = Asin 2uzo — Cs, 
which is equivalent to (C2). The result regarding the “impedance circle”’ 
follows as in Section 8. 








VEILING GLARE IN AERIAL PHOTOGRAPHY! 
By K. M. Barrp? 


Abstract 


In this paper veiling glare is defined and its general aspects in aerial photo- 
graphy are discussed. Experimental results showing the effect of glare on the 
resolving power of a standard aerial film are presented in the form of graphs. 
A method of measurement is described and typical amounts encountered in 
widely used aerial cameras are given. From the results it is shown that veiling 
glare commonly causes a loss of photographic speed of as much as 100% where 
resolving power is the primary consideration. Recommendations are made for 
minimizing flare in cameras. 


Introduction 


Veiling glare is that scattered light which causes reduction of image contrast 
in an optical system. It may be conveniently expressed mathematically in 
different ways, depending on the particular aspect studied. For example, 
L. A. Jones (3) has used a “‘flare factor’’, defined as the ratio of the luminance 
scale of the scene to the luminance scale of the image. This factor, repre- 
senting the compression of the luminance range caused by flare, was useful 
in studying tone reproduction. On the other hand, in the present paper, 
which is primarily concerned with resolving power, it has been found more 
convenient to use ‘“‘per cent veiling glare’, defined as the ratio of the amount 


of flare light in the image plane to the amount of pure image forming light, 
the ratio expressed as a per cent. 


According to the definition given above, any light scattered outside the 
geometrical image, even that caused by diffraction or residual aberrations of 
the lens, might be called veiling glare; although the term is generally used only 
to describe light scattered over a relatively large area. The latter distinction 
is arbitrary and there is no sharp division either with respect to areas involved 
or effects on resolving power. Indeed, it is probable that the results given 
below could be used to predict in lines per millimeter the resolving power of a 
lens used in combination with Aero Super XX, if the energy distribution in the 
image formed by the lens were determined by use of the same film. Such a 
method is that used by Jones and Wolfe (4). The dependence of the amount 
of veiling glare found to be present on the area considered is discussed in the 
section on measurement of veiling glare. 

In an exact study of the effect of veiling glare on the resolving power of film, 
the cause as well as the amount must be taken into account. This will be 
clear from a consideration of such typical causes as light scattered from the 
walls of the camera, reflections between lens surfaces, haze, etc. In the first 
example quoted, light that would not otherwise reach the image plane has been 
thrown on to the film, thereby increasing the illumination. In the second 


1 Manuscript received February 25, 1949, 
Contribution from the Division of Physics, Nationai Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 1952. 
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example, on the other hand, the veiling glare arises from the redistribution of 
what was originally image forming light, and the illumination of the film is not 
increased in general, although the shadow areas gain light at the expense of 
the highlights. Thus there are two general types of veiling glare, that which 
increases the amount of light reaching the film and that which does not. For 
convenience these will be designated ‘‘Type |” and “Type II’’, respectively. 
Haze does not clearly fall into one or other of these types but will be shown 
later to be more conveniently considered of Type II. The difference in 
effect of these two types will be considered in the next section. 

Of the various causes of veiling glare those due to the camera or lens are 
usually easy to classify as one or the other of the above two types. In the 
case of haze the classification is a little more difficult since in considering the 
difference between conditions when haze does and does not exist the effect 
of decreased ground illumination caused by haze must be taken into account. 

If it is assumed that an exposure meter is used as a guide to the camera 
aperture and shutter settings, it seems logical to consider the haze as veiling 
glare of Type II, i.e., it reduces contrast but does not greatly affect the average 
image luminance, since the meter reads the luminance of ground plus haze. 
Haze could be measured as a per cent veiling glare by determining the loss 
of contrast of objects that were of known contrasts on the ground. In this 
paper only measurements of veiling glare caused in the camera or lens are 
described. 

Effect of Veiling Glare 

The effects of the two types of veiling glare are indicated in Figs. 1 and 2 

as a modification of the characteristic curve of the film. The characteristics 
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Fic. 1. Characteristic curve of Super XX Aero Film showing modification caused by 20% 
veiling glare in the case where glare consists entirely of added light. The range of typical 
exposures and the total exposure range are taken from a British report (5). 
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of the film itself are unchanged of course but the reproduction of tones under 
the given conditions of veiling glare is equivalent to that by a film having the 
modified characteristic curve under conditions of no glare. For the following 
2.4 
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Fic. 2. Same as Fig. 1 except that the veiling glare arises from what was originally image 
forming light. 


discussion it is convenient to define ‘‘object exposure” as the exposure that 
would result if the causes of veiling glare were absent, as distinguished from the 
“image exposure’, which the film receives in the image plane. The corre- 
sponding terms ‘‘log object luminance ratio’ and ‘‘log image luminance ratio” 
will also be used. 

The broken lines of Figs. 1 and 2 are graphs of the densities on an Aero 
Super XX film, plotted as ordinates, that result from the log object exposures 
that are plotted on the abscissa. The solid lines of both figures are charac- 
teristic curves of Aero Super XX film and therefore identical. The broken 
line of Fig. 1 was obtained by finding the density that resulted from the object 
exposure plus 20% of the typical object exposure in aerial photography. This 
density was then plotted against the log object exposure. The densities on 
this curve are of course always greater than the Super XX densities owing to 
the added light, but the two curves approach each other asymptotically at 
large exposures. 

The broken line of Fig. 2, which represents veiling glare of Type II, was ob- 
tained by finding the density that resulted from the object exposure less 1/6 
(i.e., 5/6 of the object exposure) plus 1/6 of the typical object exposure. It is 
clear that this corresponds to 20% veiling glare as it is defined above. This 
curve has the same shape as that of Type I but differs in that it crosses the 
Super XX curve at the typical log object exposure, namely —1.4 log meter 
candle seconds. 








BAIRD: VEILING GLARE IN AERIAL PHOTOGRAPHY 133 


It should be understood that the log exposure values along the abscissa do 
not represent for the modified curves the actual log exposure on the film 
but rather the exposure that would result under conditions of no glare, i.e., 
they are log object exposures. The log image exposure corresponding to a 
given point on the modified curve is given by considering the point on the solid 
curve representing the same density as the given point. 

It is pointed out that using 20% of the exposure at the middle of the typical 
exposure range as the added light is equivalent to assuming that the average 
exposure is the same as the typical exposure. This does not exactly represent 
conditions in an aerial photograph but it is felt that the inexactness does not 
affect the conclusions drawn except when there are unusually bright sources in 
the field of view. 

The range of typical exposures and the total range of exposures shown in 
the figures were taken from a British report based on statistical findings in 
reconnaissance photography over Europe during the war (5). 

The figures show that the two types of glare have different effects. In the 
case shown in Fig. 2, at the region of typical exposures, the glare reduces the 
optical contrast without affecting the exposure, with the result that the photo- 
graphic contrast is reduced. In Fig. 1, on the other hand, the exposure is 
increased, resulting in the use of a different part of the Super XX curve. 
Since the slope of the curve is increasing, the loss of optical contrast tends to 
be compensated photographically in the lower part of the exposure range. 
It is noteworthy that in Fig. 1 the contrast on the film is actually increased 
at the toe of the curve although obviously the contrast of the optical image is 
always reduced by veiling glare. This is caused by the rapid change in slope 
at the toe of the true characteristic curve of the film. The effects on resolving 
power will be discussed in detail in a later section. 


Determination of Effect on Resolving Power 


The effect of various amounts of veiling glare on the resolving power of a 
camera-film combination was found by calculating the resulting change in 
contrast and making use of experimental curves relating resolving power to 
log exposure and contrast. A set of such curves for Super XX Aero Film 
processed in D-19) to a gamma of 1.34 is shown in Fig. 3. 

These curves were obtained by printing sets of resolution charts on strips 
of film, the charts being illuminated through a calibrated optical step wedge. 
One of the films so obtained is reproduced in Fig. 4. Resolving power readings 
from such a film, together with the calibration of the wedge, gave one of the 
curves. (Each curve in Fig. 3 represents the average of several films.) The 
. different curves were obtained by using charts of different contrasts. Proces- 
sing was done with the accurate control that is usual in sensitometry. From 
the H and D curve, which also has been plotted in Fig. 3, it is seen that the 
range of typical densities recommended in the British paper (5), quoted above, 
lies in the region of maximum resolving power for a log luminance ratio of 0. 2. 
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This value is given as the typical object contrast by the same report. In 
these curves and throughout this paper, contrast is expressed as the log of 
the luminance ratio of the two target areas. The log exposure is given for 
the brighter target area so that when the log exposure is such that the density 
resulting on the film is zero, the resolving power is also zero. 


165 















































DENSITY 








« 
wi 
é 
Oo 
z 
> 
3 
° 
” 
w 
c 
3 
S 














— AVERAGE ¥ 01.34} —7 z= 














1.5 1.0 . ° 
NEGATIVE LOG. EXPOSURE 


Fic. 3. Curves showing the variation of resolving power with log exposure for various log 
luminance ratios of the test target for Super X X Aero Film developed in D-19b. The densities 
along the right-hand side of the graph refer only to the Hand D curve. The log resolving power 
values refer to the other curves. 


The solid curves in Figs. 5, 6, and 7 were obtained from Fig. 3 by inter- 
polation and have been plotted on a larger scale for convenience. These 
graphs show three variables, giving the resolving power for any log exposure or 
contrast in the range covered. Thus by calculating the change in log exposure 
and contrast caused by a given amount of veiling glare, the change in resolving 
power can be obtained. This has been done for values of 10, 20, and 40% 
veiling glare of Type I in Figs. 5, 6, and 7 respectively. The change in 
resolving power is indicated by the arrows and by the broken curves, which 





Fic. 4. A reproduction of one of the films used to obtain the curves in Fig. 3. 


Fic. 9. Shows a picture taken in the tunnel in Fig. 8 for the determination of per cent 
veiling glare in an aerial camera. Note the grayness of the image of the black hole. 
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were obtained by calculating the modified log exposure and modified contrast 
for the given amounts of veiling glare. The modified log exposure was cal- 
culated by adding to the corresponding original exposure value a certain per 
cent of the typical exposure. By doing this for a given log exposure value 
and also for that value decreased by the log luminance ratio, two modified 
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Fic. 5. Curves showing the effect of 10% veiling glare on the resolving power of Super XX 
Aero Film. The solid curves, obtained from Fig. 3, show the variation of resolving power with 
log exposure when there is no glare. The arrows show the change in position on the graph 
caused by 10% veiling glare. The broken curves show the modification of the film characteristic 
that results. 


log exposure values were obtained, the difference of which gave the modified 
log luminance ratio contrast. 

From the modified log exposure and contrast the resolving power is given 
by the appropriate point on the graphs. Corresponding points on curves 
without veiling glare (solid) and curves with veiling glare (broken) are repre- 
sented by the tail and the head of an arrow respectively. Hence the change 
in resolving power due to a certain per cent veiling glare may be found by 
reading the head and tail of an arrow. For example log object exposure 
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= —1.4 and log object luminance ratio = 0.2 gives a log resolving power 
of 1.4 if there is no veiling glare. In the presence of 20% veiling glare of 
Type I this is changed to 1.35. 
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Fic. 6. Sameas Fig. 5 for 20% veiling glare. 
5 J /€ 53 


The effect of Type II veiling glare can also be obtained from the graphs. 
It can be shown that in this case the modified contrast is the same as for Type | 
but the modified log exposure is less than for the latter by an amount log 
(1 + V/100), where V signifies the per cent veiling glare. Thus for the point 
considered in the preceding example the modified log exposure is the same as 
the original (as would be expected at the middle of the typical exposure range) 
and the modified resolving power is therefore again about 1.35 (given by a 
point slightly above the broken curve). 

The graphs show that Type II veiling glare caused greater loss of resolving 
power at low exposures but less loss at high exposures than Type I. However, 
as the above example shows, the difference is not great for the log exposure 
range of —1.3 to —1.5, i.e., in the range of typical exposures in aerial photo- 
graphy. Infact the difference is not significant for the arguments in this 
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paper, and therefore the two types will not be considered separately hereafter. 
In any case, in practice both types exist simultaneously. 
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Fig. 7. Sameas Fig. 5 for 40% veiling glare. 









Discussion of Results 






Considering a log luminance ratio of 0.20, Fig. 5 shows that 10% veiling 
glare changes the log resolving power by about —0.03, i.e., a reduction of 
about 7%, in the typical exposure range. This amount is of the same order 
as the variation in lens-film resolving power from lens to lens in good quality 
survey lenses (2) and is therefore not serious. On the other hand 20 and 40% 
cause losses of about 12 and 20° respectively. The latter amounts are great 
enough to be undesirable but frequently occur in practice, as will be seen 








in the next section. 





Instead of considering the effect of veiling glare as a loss of resolving power 
it may be thought to constitute a loss of camera speed. Thus if film speed 
is determined using the resolving power criterion suggested by L. E. Howlett 
(1), then the minimum acceptable exposure (i.e., that giving 90% of the 
maximum resolving power when the log luminance ratio is 0.20) is —1.70 
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for the film being considered, in the absence of veiling glare (see Fig. 6). 
To realize the same log resolving power (viz., 1.35) in the presence of 20% 
veiling glare on the other hand, the log exposure must be —1.35. This is 
an increase in log exposure of 0.35 or an increase in the exposure of over 100%. 
As a result a camera equipped with an f/5.5 lens having 20% veiling glare 
would be effectively no faster in cases where resolving power is the prime 
consideration than one with an f/8 lens free from glare, assuming the same 
quality in other respects. While this admittedly may exaggerate the situation 
(for example if the minimum criterion had been 95% of the maximum, the 
loss in speed would have been infinite) it is felt nevertheless that the effect 
is more serious than is generally realized. To state the effect generally: for 
the useful range of exposures any resolving power obtainable in the presence 
of veiling glare can be achieved with less exposure if the glare is reduced. 
This effective increase of camera speed is more significant than gain of trans- 
mission if lens surfaces are coated to reduce reflection. An example is 
described in the section ‘‘Veiling Glare in Some Existing Cameras’’ that 
follows. Corresponding changes of photographic speed as determined by 
density criteria can be seen by considering Figs. 1 and 2. 

It is interesting that for very low exposures 10 or 20% veiling glare increases 
the resolving power. It is this effect that is sometimes used in spectrographic 
work to detect very faint lines by deliberately fogging the photographic plate. 
However, over the range of exposures used in aerial photography the glare 
causes a loss of resolution. 

Since the change of resolving power has been determined for the film alone, 
it might be argued that the results would not be valid for a lens—film com- 
bination. However it is felt that the conclusions would be the same in the 
latter case because the lens-film resolving power does not generally differ 
greatly from the resolving power of the film alone when fast relatively coarse 
grained film is used, as is the case in aerial photography. 


Measurement of Veiling Glare 


The amounts of veiling glare caused by inherent defects in camera-lens 
combinations have been determined for several widely used aerial cameras. 
The method of measurement was to photograph a field of view, which was of 
uniform luminance except for a black spot subtending an angle of about one 
degree, with the camera under test. A specially constructed tunnel shown 
diagrammatically in Fig. 8 was used to provide this field. 


In the film plane of the camera a small calibrated optical step wedge was 
placed so that part of the bright field became exposed through it. An example 
of the pictures that resulted is reproduced in Fig. 9. After processing the film 
carefully, using brush development, suitable temperature control, etc., densities 
were read for the image of the steps of the printed wedge, the black spot, and 
the field. The characteristic curve of the film was plotted from the wedge 
densities and .the log exposure difference between the images of the black spot 
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and the background obtained from the curve as indicated in Fig. 10. The 


, j ; ‘ Glare : P 
antilog of the value obtained is the ratio: ~———-—_——.._ -_ If this ratio 
Glare + Image light 







; ss ats R 
is called R then the per cent veiling glare is given by: oe x 100. 
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Fic. 8. Diagrammatic view of the tunnel used to measure veiling glare in aerial cameras. 
The sections of the tunnel containing the lights are painted white; the rest is black. 
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Fic. 10. Characteristic curve of the film shown in Fig. 9 plotted from the wedge included 
in the picture. The dotted lines indicate how the log exposure difference between the image of 
the black hole and the background is determined. 










The value for ‘Glare + Image light’? was measured in the immediate 
neighborhood of the point in the field where the veiling glare was being deter- 
mined. This is in general different from a value averaged over the whole 
field of view because of vignetting, the cosine effect, variation over the field of 
the amount of flare, and so on. The choice of the former value is obviously 
the logical one in view of the method of determining effects on resolution 







discussed in the previous sections. 





Recalling the discussion in the second paragraph of the introduction it is 
seen that, because the black spot used in the measurement described above 
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subtends about a degree, only light that is thrown at least a half degree from 
the edge of the geometrical image is included in the value obtained for the 
glare light. This is an arbitrary amount, and no attempt is made to justify it 
in this paper except to say that veiling glare has been used to describe only 
that light which is scattered over a ‘‘relatively large’ area. In the case where 
secondary images are formed by reflection between lens surfaces it is clear 
that the scattered light may involve any area, as the secondary image may be 
very much out of focus or in focus and coincident with the primary image. 

It is apparent that the method described above zives the per cent veiling 
glare that is effective when the camera is used in service because the color 
sensitivity, diffuseness of pickup of energy, etc., is identical in the test and in 
service use, provided the same kind of film is used. Conditions have necessarily 
been idealized in that a uniformly bright field of view has been used but this 
does not lead to serious error except in special circumstances, such as cases 
when the sun is in or nearly in the field of view, which cases must be con- 
sidered separately. 

The usual objections to photographic photometry such as reciprocity law 
failure, variation in processing, etc., are obviously not valid here because of 
the method used. Reproducibility of measurements has been found consistent 
with the,practical considerations involved. For example, three separate tests 
on the axis of Camera D in Table I gave the values: 20, 25, 24% respectively. 


Veiling Glare in Some Existing Cameras 


The glare found in several widely used aerial cameras is shown in Table I. 
A and B, a short focal length reconnaissance camera and a survey camera 
respectively, are representative of cameras having fair correction. | The cones 
are painted black inside and are large enough so that there is little tendency 
for light to be reflected specularly on to the film. Neither camera has coated 


TABLE I 


Veiling glare on axis, | Veiling glare at edge 
% of field, % 


Noe wm 
NWP Ub Ww 





The identification of the cameras mentioned in the preceding table is as follows: 
A isan F24 reconnaissance camera. 
Bisa Fairchild K 17b with a 6 in. Bausch & Lomb Metrogon lens. 
Cand C’ represent a Williamson O.S.C. with a 6 in. Ross W.A. Survey lens. 
D_ isa Fairchild F224 Cartographic camera with a Metrogon lens. 
E isan F52 camera with a 36 in. F/6.3 N.R.C. Booth Telephoto lens. 





BAIRD: VEILING GLARE IN AERIAL PHOTOGRAPHY 141 


lenses, but the curvatures of the lens surfaces of Camera A are not very steep. 
Camera B has steeply curved lenses, and coating the lenses would likely 
reduce the glare to a value close to that of Camera A. The results given in 
the preceding section indicate that these cameras have lost about 10% of their 
potential resolving power because of veiling glare. This amount is not serious 
although it would be worthwhile to reduce it, especially in Camera B. 

Camera C is a recently designed aerial survey camera that is well baffled. 
The unusually high veiling glare was found to be caused by reflections from 
the lens surfaces, which are very steeply curved, with the result that reflections 
at large angles of incidence cause light to be thrown back into the camera. 
C’ is the same camera after the lens was coated with antireflection films. The 
glare has been reduced to an unobjectionable amount. 

Another recently designed survey camera with not much more than normal 
glare under normal conditions is Camera D. However, apparently in the 
interests of compactness, the cone was designed to fit the outside rays very 
closely, with a disregard for veiling glare, and no field limiting baffle was 
provided outside the lens. The result is slightly increased glare (especially 
at the edge of the field) under normal conditions, viz., when the field of view 
is of nearly uniform luminance. When very bright sources such as clouds or 
lakes with the sun reflected from them are imaged on the inside of the cone, 
serious glare results. From measurements on the camera and on aerial nega- 
tives showing lakes it was found that glare of over 35% could be expected 
frequently. 

Camera E is a well designed long-focal-length reconnaissance type with 
relatively flat lens curves. Fig. 11 shows typical variation of glare with 


angle in the field. 


% VEILING GLARE 
= - yt nN 
uw ° w ° Ww 


° 


° ' 2 3 4 5 
DISTANCE FROM CENTER OF FIELD (INCHES) 


Fic. 11. Curve showing typical variation of veiling glare across the field of view of an aerial 
camera. 


Conclusion 


From the foregoing results it is seen that in an aerial camera in which 
reasonable care has been taken in design, veiling glare will be of the order of 
10 to 15%, corresponding to a reduction in resolving power in the typical ex- 
posure range of about 10%. Whereas this amount is undesirable it is not 
serious and it is generally not practical to reduce it. 
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It is also evident however that lack of care in design might result in veiling 
glare of the order of 40%. This would cause a loss of resolving power of over 
20% or an equivalent serious loss of speed. 

To keep veiling glare at a minimum it is recommended that the lens surfaces 
be coated to’reduce reflection if there are many elements or if the surfaces 
are steeply curved, as is usually the case with wide angle lenses. There should 
be a field limiting cone outside the lens. The lens tunnel should be large, 
well baffled, and blackened on the inside. 

Whereas the experimental results given above are for aerial cameras and 
one photographic emulsion, the general conclusions apply to photography in 
general. 
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RELATIVE YIELDS OF THE D-D REACTIONS AT LOW 
ENERGIES! 


By T. P. PEPPER 


Abstract 


By the simultaneous detection of the p, T, and He’ particles in a thin window 
proportional counter the concurrent nuclear reactions D(d,p)T and D(d,n)He* 
have been investigated throughout the deuteron energy range 26 to 63 kev. The 
ratio of the yield of the reaction D(d,p)T to that of the reaction D(d,n)He’ was 
found to be 1.15 + 0.15 throughout this energy range. There is some slight 
evidence for an increase in this ratio towards lower bombarding energies. 


Introduction 


The bombardment of deuterium with deuterons produces the concurrent 
nuclear reactions , 
D(d,p)T Q = + 4 Mev. 
and 


D(d,n) He’ Q 


ll 


+ 3.3 Mev. 


which will be referred to collectively as the D-D reactions. Various aspects 
of these reactions have been reported in the literature. When the present 
work was done early in 1948, however, the literature contained no reference 
to the simultaneous measurement of the D—D reactions at Jow energies. In 
addition it appeared that no study of the yields had been made using similar 
methods of monitoring the two reactions, as, for example, the detection of 
the charged products from each reaction as herein described. The present 
work was undertaken to show that the yields from both reactions could readily 
be measured simultaneously and under identical conditions. This method 
could therefore provide a simple means of obtaining the ratio of the yields at 
low energies. 

Reasons for the present experiment are suggested below. Measurement 
of the reactions D(d,n)He*® by means of the He’® particles would provide a 
useful standard of neutron emission. Obtaining the relative yields of the 
two D-D reactions would permit a neutron standard to be established from 
measurements of the yield of the more readily detectable protons and tritons. 
Knowledge of the ratio of the yields of the D-D reactions as a function of 
bombarding deuteron energy might prove useful in estimating the relative 
importance of the Coulomb and centrifugal barriers in the formation of the 
reaction products. 


1 Manuscript received March 10, 1949. 


Contribution from the Nuclear Physics Branch, Chalk River Laboratory, Division of 
Atomic Energy Research of the National Research Council of Canada. Issued as N.R.C. No. 1956. 
This work was part of a thesis submitted to McGill University in partial fulfillment of the require- 
ments for the degree of Doctor of Philosophy, October 1948. 
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Experimental Method 


A small 50 kv. accelerator, designed and built by D. Roaf, employed a 
conventional 25 kv. voltage-doubler circuit. The ion source was a radio- 
frequency positive ion source designed and built by A. J. Bayly and A. G. 
Ward (1). The ion beam was employed without resolution. It consisted, as 
described by Ward and Bayly, of about 200 wa. mixture of masses 1 to 6, 
the odd integral masses being due to hydrogen impurities in the deuterium. 
Approximately 50% of the total beam was atomic heavy hydrogen: the relative 
intensities of D+, DD+, and DDD? were 6, 2, and 1: the relative energies were 
E, E/2,and E/3. The yields at energy E were one to two orders of magnitude 
greater than those at E/2 or E/3. To a good approximation, therefore, the 
observed yields were due to the higher energy, more intense atomic deuteron 
beam. It is the energy of this beam that is quoted throughout. For extract- 
ing and focusing the ions two separate 0 to 10 kv. power supplies were available. 
Thus an ion had a maximum energy of 70 kev. The voltage was measured 
directly by observing the current through a 50 megohm chain of carbon 
resistors. Since the present method of determining the relative yields of the 
D-D reactions did not require precise determination of this voltage, the resistor 
chain was not calibrated for voltage or thermal variations. 
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Fic. 1. Cross sectional view of a water-cooled copper target for condensing the deuterons and 
proportional counter for detecting protons, tritons, and He particles. 


Fig. 1 represents a cross-sectional view of target and cylindrical counter. 
The deuterium target was a condensation target, i.e., a water-cooled copper 
target onto which the impinging ions condensed, forming a deuterium layer 
that was then bombarded by succeeding deuterons in the beam. The result- 
ant particles from the D-D reaction, i.e., 3 Mev. protons, 1 Mev. tritons, 2.4 
Mev. neutrons, and 0.8 Mev. He® particles, penetrated a thin Formvar 
window and entered the sensitive volume of the proportional counter in a 
cone of semiangle’20°. This degree of collimation was effected by an array 
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of 0.040 in. holes drilled in a 0.1 in. steel plate, which supported the Formvar 
window. During the experiment, windows of 0.25 and 0.12 mgm. per cm.? 
were used. A thin layer of gold was sputtered onto the Formvar windows to 
render them conducting and to reduce the light entering the counter from 
the deuteron beam and target. The minimum path length in the counter 
for the charged particles was 3cm. For most of the runs, methane was used 
as counter filling. The final series of runs was made, however, using -octane 
vapor at room temperature (7). The counter was operated with a gas ampli- 
fication of approximately 100. The ionization pulses were further amplified, 
sorted according to size, and counted. 


Experimental Results 


The experimental procedure was arranged, first, to identify the three 
groups of pulses observed on the pulse analyzer (8), and, second, to obtain the 
relative yields of the D-D reactions. The results are most readily presented 
by referring to representative pulse distributions observed. 
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Fic. 2. Proton, triton, and He* pulses from the D-D reactions. Methane-filled propor- 
tional counter at pressure of 30 cm. of mercury. Formvar window 0.25 mgm. per cm.? 
Energy = 39 kev. 


Figs. 2, 3, and 4 show the three groups of pulses labelled as subsequently 
identified. In Fig. 2 (counter pressure, 30 cm. of mercury) the He* peak is 
broad and spreads across the sharper, less energetic proton peak, while the 
triton peak falls entirely in Channel No. 30, the last channel, of the pulse 
analyzer. This channel counted all pulses above a certain size. By reducing 
the pressure in the counter the triton peak emerges as shown in Fig. 3 (counter 
pressure, 18 cm. of mercury) to become a well resolved peak while the He*® peak 
is somewhat sharpened and the proton peak has receded into the noise back- 
ground of Channel No. 1. By stil! further reducing the pressure in the counter 
the relative positions of the triton and the He* peaks were interchanged, as 
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expected (Fig. 4). At counter pressures of 1.9 cm. of mercury the He* peak 
became sharp and similar to the proton or triton peak. 
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Fic. 3. He®and triton pulses fromthe D-D reactions. Methane-filled proportional counter 
at pressure of 18 cm. of mercury. Formvar window 0.25 mgm. per cm.2 Energy = 60 kev. 
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Fic. 4. Triton and He’ pulses from the D-D reactions. Methane-filled proportional counter 
at pressure of 3.5 cm. of mercury. Formvar window 0.12 mgm. percm.2 Energy = 43 kev. 


The ratio of the yields of the D-D reactions was obtained in three ways 
from the resolved proton, triton, and He* peaks. These are best illustrated 
with reference to Figs. 2 and 4. 

Method A 

For the counter filled with methane to a pressure of 30 cm. of mercury 
(Fig. 2) the triton pulses fell in Channel No. 30. By counting these pulses 
and the pulses associated with the combined proton and He? peaks, the ratio 


T 


of the number of counts —————————=. can be obtained. This is equivalent 
(i (p + He’) — 1 ’ 
to the ratio T/He’. 
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Method B 
Again from Fig. 2 the ratio p/He* can be obtained by direct interpolation 
of the He* peak beneath the p peak. 


Method C 

For the counter filled with methane to a pressure of 3.5 cm. of mercury 
the T and He? peaks were resolved as illustrated in Fig. 4. Upon an arbitrary 
assignment of those counts lying between the main portions of the T and He* 
peaks a straightforward deduction of the ratio T/He* is obtained. Two 
extreme cases were considered: the one apportioned the ambiguous counts 
so as to maximize the ratio T/He’, the other to minimize it. The largest 
spread encountered in all the cases considered was 14%. The consistent 
cleanness of the T peak and the characteristic straggling of the He* peak 
lead to an arbitrary weighting of 2:1 in favor of assigning the doubtful 
counts to the He? peak. 

Ratios obtained by Methods A, B, and C are each treated as equivalent to 
the ratio of the yields of the reactions D(d,p)T and D(d,n)He*. Had the 
peaks been completely resolved from one another and any spurius background, 
all three methods would have given the same results. The three methods are 
presented and employed therefore to indicate the amount of uncertainty 
introduced by the imperfect resolution. Averaging of the results removes in 
large part the errors associated with any one of the methods outlined above. 

Each method was applied at various deuteron energies from 26 to 63 kev. 
The ratios so obtained are presented in Table I. 


TABLE I 


MEASUREMENT OF YIELDS OF D(d,p)T AND D(d,n)He® 


d-energy Counter filling, * Ratio of 
kev. in cm. Hg yields 


Method A 
Ratio calculated 27 : 1.07 
= 1.06 


(p + He’) — T 1.10 


3 ‘ 0.99 
Method B ‘ 1.07 


Ratio calculated 45 ‘ ‘81 
p 4 1.12 


3 1.12 
a 1.12 





1.35 
1.19 


Method C 3.5 1.20 
Ratio calculated Si 


He? 
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Discussion of Results 


Method C was applied to runs made using a window 0.12 mgm. per cm.?, 
i.e., half the thickness of the one providing the data for Methods A and B. This 
was permitted by the lower pressure requirements of Method C. That the 
ratios obtained by Method C are some 8% higher than those obtained by 
Methods A and B may indicate selective absorption of the He* particles by 
the window. The ratios 1.31 and 1.28 deduced using Method C at bombard- 
ing energies 26 and 61 kev. were obtained using m-octane vapor at room tem- 
perature. These ratios are in satisfactory agreement with those obtained 
using methane at about six times the pressure. 

The possibility of selective absorption of the charged products by small 
irregularities in the surface of the target was considered but thought to be 
unlikely, as indicated by the following argument. To the deuteron beam 
the target depth was just that for the energy of the deuteron beam. This 
depth / was never greater than 63 kev. The target made an angle of approxi- 
mately 30° with the deuteron beam. The particles entering the counter 
therefore penetrated a maximum target of / tan (30 + 20)°. For a perfectly 
smooth target only the shortest range particles (He* of 4 mm. air equivalent) 
need be considered. The He* particles, would, on their way to the counter, 
lose less than 100 kev. in the target. With the windows used this would 
broaden but not decrease the He* peak. On the other hand, for a target with 
small metal spurs, the emitted particles might indeed be prevented from 
entering the counter. However, if such were the case for an appreciable 
fraction of the particles, those that did enter the counter would not form a 
defined peak. The fact that the He* particles did produce such a peak, 
which with adequate reduction of the counter stopping power could be 
sharpened to resemble the p or T peak, is good evidence against serious loss 
of the He’ particles before entry into the counter. Asa precautionary measure 
the copper target was always very highly polished with Hydro-Durexsil 700 
paper and only in an up and down motion so that whatever grooves were 
formed would be in the plane containing the beam, the target, and the counter 
window. 

The possibility that the measured ratio was altered by scattering of the 
emitted particles on their way to the counter was considered. Assuming the 
scattering to be of the same order of magnitude as that of deuterons on deute- 
rium as reported by Blair et a/. (2), the fraction of p, T, He’ particles scattered 
before entry into the counter was found to be negligible. 

The D-D reaction products are known to be anisotropically distributed 
in the center of mass system of co-ordinates. Bretscher, French, and Seidl (4) 
have measured the angular distribution of protons for deuteron energies 20 to 
80 kev. They found that the angular distribution could be represented by 
the relation NV, = Ngo(1 + A cos? 8), where N, is the proton flux proceeding 
at an angle @ to the deuteron beam and Np is the flux proceeding perpendi- 
cularly to the beam. The quantity A is a function of the deuteron energy 
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and increases linearly from 0.25 at 20 kev. to 0.45 at 80 kev. The angular 
distribution of the neutrons has not been determined at these low energies. 
At higher energies, however, the angular distributions are similar for protons 
and neutrons (3). It is therefore unlikely that they differ radically at low 
energies. In the case of identical angular distribution functions for protons 
and neutrons the ratio of the fluxes at 6 = 90° can correctly be taken as the 
ratio of the total yields of the reactions D(d,p)T and D(d,n)He®. In the 
extreme case, that the neutron flux is isotropic for the energies considered, 
the error so introduced is only 6%. That the angular distribution of the 
neutrons is markedly more asymmetric than that of the protons is very 
unlikely. It is improbable therefore that the error introduced by the angular 
distribution is more than 2%. 

From Table I it is seen that the ratio of the yields of the reaction D(d,p)T 
and D(d,n)He? is 1.15 + 0.15 throughout the energy range 26 to 63 kev. 
Assuming that the results obtained using the thinner window and Method C 
are the more reliable there is some slight evidence that the ratio increases 


toward lower deuteron energies. 

Blair et al. (3) have found that throughout the energy range 1 to 3.5 Mev. 
the ratio of the cross sections of the reactions D(d,p)T and D(d,n) He’ is almost 
constant and equal to 0.9. The reaction D(d,n)He’* was detected by means 
of the He’ in a proportional counter, but owing to the high energy of the 


scattered deuterons the He’ particles and protons were not monitored simul- 
taneously but at different times using windows of different thickness. At 
the low deuteron energies of the present experiment, the scattered deuterons 
could not enter the counter, and hence the charged products could readily be 
detected together. 

Coon and co-workers (5, 6) have investigated the D-D reaction throughout 
the energy range 100 to 300 kev. and have found that the ratio of the cross 
sections of the reactions D(d,p)T and D(d,n)He?’ is 0.6 at 300 kev. and 1.0 
at 130 kev. Their results indicate by extrapolation that at still lower energies 
the ratio exceeds unity. These results were obtained by detecting the protons 
in an ionization chamber at angles of 20°, 45°, and 90° to the deuteron beam 
and the neutrons in a surrounding manganese sulphate solution. The relative 
cross sections therefore required the use of a neutron standard, a complication 
not met when the reactions are detected entirely by the charged particles. 
The present results are consistent with those of Coon et al., extrapolated to 


lower energies. 
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THE SPECIFIC HEAT OF NAPALM*-GASOLINE GELS! 


By G. O. LANGSTROTH? AND K. H. Hart 


Abstract 


The heat capacities of gasoline, Napalm, and various Napalm-gasoline gels 
have been measured in the range —50° C. to 50° C. by the method of mixtures, 
using an ordinary type of Richard's adiabatic calorimeter. The specific heat in 
calories per gram per Centigrade degree at a temperature T° C. within this range 
is given to within 2% by the relation 


Sr = 0.479 — 0.00054C + 0.000927, 


where C denotes the Napalm concentration in per cent by weight. Concentra- 
tions of 0% and 100% are included. These data supplement earlier results on 
thermal conductivity. 


The thermal conductivity of Napalm-—gasoline gels has been determined 
previously (2) over the temperature range —50° to 50°C. with Napalm 
concentrations from 0 to 10% by weight. Specific heat data for the same 
temperature range with Napalm concentrations from 0 to 100% are presented 


in this article. 
Apparatus and Procedure 


The determinations were made by the method of mixtures with the ordinary 
type of Richard’s adiabatic calorimeter. The cylindrical brass calorimeter 
vessel (diameter, 2 in.; height, 4 in.) was mounted on an insulating support 
inside a water-tight brass jacket (diameter, 3 in.; height, 53 in.). The jacket 


was immersed in a motor-stirred water bath of 3 gal. capacity. Samples were - 
introduced to the calorimeter vessel through a tube soldered to the top of the 
jacket and extending just above the bath level; normally the tube was closed 
with a cork. The calorimeter vessel was provided with a hinged top, which 
was opened or closed by a thread extending through the tube. The water in it 
was stirred manually by means of a small copper vane mounted on a plastic 
rod that passed through the top of the jacket. 

A copper-constantan thermocouple, with one junction in contact with the 
calorimeter vessel and the other in the bath, was connected directly to the 
terminals of a galvanometer. A difference of 0.003° C. between the calori- 
meter and bath temperatures produced a deflection of 1 mm. Changes in 
bath temperature involved in maintaining adiabatic conditions while the 
sample and calorimeter contents attained equilibrium were read to 0.001° C. 
from a Philadelphia differential thermometer. Bath temperatures were read 
to 0.01° C. from a calibrated mercury thermometer. 

Samples to be measured were placed in cylindrical containers of thin brass 
closed by screw caps with thin gaskets. The containers had a diameter of 


1 Manuscript received March 3, 1949. 
Contribution from the Physics Department, University of Alberta, Edmonton, Alta. 
2 Present address: Experimental Station of the Defence Research Board of Canada, Suffield, 
Alta. 
* Napalm is a commercial product consisting of a basic aluminum soap of naphthenic, oleic, 
and palmitic acids in the approximate ratiol: 1: 2. 
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2 cm., a height of 7 cm., a volume of 20 €c., and a mass of 18 gm. The gasket 
weighed about 0.1 gm. No leakage was detected. 

The bath used to bring the samples to the desired initial temperature before 
introduction to the calorimeter had a capacity of about 1 gal. It was similar 
to that used previously in conductivity measurements (2). The bath liquid 
was water for temperatures above, and methyl alcohol for those below, 0° C. 
Temperature constancy was maintained by the operation of a heating element 
with a toluene thermostat control working against a cooling tube containing 
acetone and carbon dioxide snow. Vertical brass tubes about 7 in. long and 
1 in. in diameter with their tops just above the liquid level served as dry 
receptacles for the sample containers. The bath temperature was measured 
to 0.01° C. by potentiometric determination of the e.m.f. of a copper—con- 
stantan thermopile of three couples in series, calibrated at the steam, ice, and 
carbon dioxide snow points. It exhibited periodic fluctuations of 0.01° C. 
above 0° C., and of 0.05° C. below 0° C., but no drift occurred over periods 
of several hours. 


Before insertion in the constant temperature bath, the sample containers 
were placed in shields designed to minimize heat losses or gains during the four 
second interval required for transfer from the bath to the calorimeter. A 
shield consisted of a brass tube just large enough to hold a container. It was 
closed permanently at the upper end and was corked at the lower. In trans- 
ferring a sample from the bath to the calorimeter, the shield was held over 
the tube in the calorimeter jacket, the cork was removed, and the sample 
container was lowered into the calorimeter by means of a thread passing 
through the top of the shield. 

A graded series of gels was prepared by mixing appropriate quantities of 
Napalm and gasoline with mechanical stirring for half an hour. The gasoline 
was supplied by the Research Council of Alberta. It had a tetraethyl lead 
content of 3.4 cc. per gal., and a specific gravity of 0.702. About 15 gm. of 
the gel to be examined was placed in a sample container and brought to the 
desired initial temperature by insertion for at least one hour in the constant 
temperature bath. The temperature attained by it was taken to be equal 
to the mean bath temperature over the latter part of this period. The 
temperature of the weighed mass of distilled water in the calorimeter vessel 
and also that of the calorimeter bath was so adjusted that the value attained 
after adding the sample was as nearly as possible 21.00° C. On transfer of 
the sample to the calorimeter, adiabatic conditions were achieved as nearly 
as possible by additions of warm or cold water to the calorimeter bath. On 
attainment of thermal equilibrium the temperature change of the bath, and 
so of the water in the calorimeter vessel, was determined from the initial and 
final readings of the differential thermometer. The final temperature was 
read to 0.01° C. from the bath thermometer. The heat capacity of the sample 
between its initial and final temperature was calculated in the usual way, the 
heat capacities of the calorimeter vessel, stirrer, and sample container being 
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taken into account. The heat capacity between the initial temperature and 
21.00° C. was calculated from this value by application of a small correction. 

From 100 to 140 gm. of water was used in the calorimeter vessel; its tem- 
perature never exceeded 24° C. nor fell below 18°C. Temperature changes 
caused by introduction of the sample were from 2° to 4° C. The specific heat 
of water was taken as 1.000 cal. per gm. per Centigrade degree. In view of 
the small departures from 21° C., use of this value introduced an error of 
less than 0.1%. 

The heat capacity of the calorimeter vessel and stirrer was determined in 
tests in which 5 to 8 gm. of distilled water at 0° C. was added to 20 to 30 gm. 
of water in the vessel, with a resulting temperature drop of about 3° C. The 
mean value for the 21°C. region was 16.5 + 0.2 cal. per °C. The heat 
capacity of each sample container was calculated from the masses and specific 
heats of brass and the gasket material. On the average these contributed 
respectively about 0.05 and 1.6 cal. per °C. to the heat capacity of the loaded 
container (9 to 10 cal. per °C.). The specific heat of the gasket was deter- 
mined from measurements on sample containers filled with disks of the 
material. It was found to be 0.52 cal. per gm. per Centigrade degree. Since 
the gasket accounted for only about 0.5% of the heat capacity of the loaded 
containers, variations with temperature were not measured. The specific 
heat of commercial brass was determined from heat capacity measurements 

‘on a 150 gm. cylinder by the method described above. The data are given 
in Table I. The heat capacities Hx between the initial temperature # and 
21.00° C. are described to better than 1% by the relation 


Hn = — 1.859 + 0.08854 + 0.00001962 (1). 
Consequently the specific heat Sr at T° C. is given in calories per gram per 
Centigrade degree by Sp = 0.0885 + 0.0000387 (2) 
over the temperature range of Table I. 


TABLE I 


OBSERVED HEAT CAPACITIES IN CALORIES PER GRAM BETWEEN @ AND 21.00° C. FOR A SAMPLE 
OF COMMERCIAL BRASS ROD, AND THE DEVIATIONS FROM THEM OF THE VALUES GIVEN 
BY EQUATION (1) 


8 ¢(E,) Cal. /gm. % Dev. 





.63 ; 27 —3.74 
49.64 .59 F 85 —4.49 
40.18 ota : 39.76 —5.33 
33.40 11 ; | <oS —6.19 
0.00 85 en | Kae —6.32 


In preliminary tests the bath temperature was maintained equal to that of 
the water in the calorimeter over periods of more than half an hour without 
detectable change in the reading of the differential thermometer. It follows 
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that heat exchange with the exterior and heat loss by evaporation were 
negligible over the period required for making a determination. Other tests 
showed that 100 rather vigorous strokes of the stirrer caused a temperature 
rise of only 0.005° C. in the water in the calorimeter. Since in making a 
determination only three or four gentle strokes were used on about three 
occasions during the course of heat exchange, the stirring was not sufficient 
to cause a detectable rise in temperature. As a check on over-all reliability, 
the average specific heat of copper was determined for the ranges 21° to 50° C. 
and 21° to 58°C. The mean of six determinations was 0.0932 cal. per gm. 
per Centigrade degree, in agreement with a value of 0.0930 cal. per gm. per 
Centigrade degree calculated from the results of Bronson et al. (1). 















Results of Heat Capacity Measurements 

Measured heat capacities between various initial temperatures @ and 
21.00° C. are given in Table II for gasoline, Napalm, and gasoline-Napalm 
gels of four different concentrations. 








TABLE II 


OBSERVED HEAT CAPACITIES IN CALORIES PER GRAM BETWEEN 9 AND 21.00° C. FOR GASOLINE, 
NAPALM, AND VARIOUS GASOLINE-NAPALM GELS, AND THE DEVIATIONS FROM THEM OF 
THE VALUES GIVEN BY EQUATION (4). PERCENTAGE NAPALM CONCENTRATIONS BY 

WEIGHT ARE GIVEN AT THE HEAD OF EACH PAIR OF COLUMNS. 












12. 





5 












































14.87 
49.64 14.79 | -—0.8 - - — — 
44.61 12.05 |—0.4 11.89 0.3 11.96 | -—0.2 11.91 0.3 11.74 0.9 10.72 0.7 
40.18 9.84 |—-1.3 _— _— _ — - - - : 
39.81 = — —- _ _— 9.43 0.3 _ _ 8.47 4.2 
39.73 9.56 |—0.9 9.38 0.4 9.20 ace —_— —_ 9.40 |-—0.4 — _ 
33.42 6.24] 0.3 a = —_ — — a — a 
30.52 22: } 2A 4.62 2.8 4.63 1.2 4.63 2.4 4.79 |-1.7 4.38 |-1.1 
0.00 |—10.47 |—1.9 |—10.07 1.4 |-—10.18 0.1 |}—10.10 0.5 |— 9.95 1.7 |— 9.03 0.4 
—11.17 |—15.66 | —0.6 |—15.41 0.4 |-—15.35 0.5 |-—15.47 |—0.6 | —15.30 0.2 |—13.88 | -—0.9 
—21.27 |—20.22 0.2 | —20.01 0.6 | —20.19 |—0.4 | —20.03 0.0 | —19.99 |—0.2 |—18.00 | —0.5 
—29.85 —_— — —_ _ — —_— —23.89 |-—0.1 _ _ —21.54 |-—0.8 
—31.26 |—24.85 |-—0.2 | —24.49 0.7 |—24.65 |—0.2 _ _ —24.40 0.2 — _ 
—39.76 |—28.59 |-0.1 | —28.32 0.2 | —28.48 0.5 | —28.28 |—0.1 | —28.23 |—0.2 | —25.25 0.0 
—49.53 | —32.85 0.1 | —32.38 0.9 | —32.53 0.2 | —32.45 0.1 | —32.36 0.1 |—29.08 |—0.3 
—50.77 | —33.62 ;—0.7 | —33.06 0.5 |—33.23 |—0.3 | —33.12 |—0.3 | —33.11 |-—0.5 |—29.50 |-0.1 
Average ee ue — 0.8 —_ 0.7 —_ 0.5 _ 0.6 _— 0.7 











































Discussion 


An examination was made of two types of empirical relation to describe the 
heat capacity data of Table II. The one finally adopted was the power series 


Hn = — P+ Q0+ R@, (3) 













LANGSTROTH AND HART: SPECIFIC HEAT OF NAPALM-GASOLINE GELS 155 


where Hs, denotes the heat capacity in calories per gram between the initial 
temperature 6 and 21.00°C., and P, Q, and R are constants. The values 
of the constants were first determined by fitting the equation by the method 
of least squares to the data independently for each Napalm concentration, 
including 0% and 100%. Plots of P, Q, and R against the Napalm concen- 
tration C then showed a linear relation between the first two and C, but no 
detectable dependence on C of the small constant R. The equations of the 
lines in the P and Q plots were substituted for P and Q in Equation (3) to 
express explicitly the dependence of Hy, on Napalm concentration. The 
value of R was taken to be the mean value for the entire data. On this basis 
Equation (3) became 


Hy = — 10.27 + 0.012C + (0.479 — 0.00054C) @ + 0.00046 62. (4) 


The values of Hy calculated from Equation (4) showed an average deviation 
of 0.7% from the 69 observed values, as mav be seen from Table II. The 
maximum deviation was 2.8%. There were three deviations greater than 
2%, and 11 between 1% and 2%. Two of the former, and four of the latter 
were associated with data for the initial temperature closest to 21.00°C., 
and for these the experimental error is expected to be greatest; the average 
deviation for the 6 = 30.52° data was 1.8%. Equation (4) satisfies the 
requirement that H3, should effectively vanish when @ = 21.00°C. Cal- 
. culated values of Hj for the concentrations in the order given in Table II 
are 0.03, —0.01, —0.02, 0.00, 0.00, and 0.06 cal. per gm. 


It follows from Equation (4) that the specific heat S7 in calories per gram 


per Centigrade degree at any temperature T° C. between 50° and —50° C. 
and any Napalm concentration C in per cent by weight, is given by 


Sr = 0.479 — 0.00054C + 0.000927. (5) 


This relation is considered to be good to about 2%. Although no data were 
obtained between C = 12.5% and C = 100%, there is no reason to expect 


4 


serious errors from interpolation for this region. 
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THE NUCLEAR MOMENTS OF P*'! 


By M. F. CRAWFORD AND J. LEVINSON? 










Abstract 


Hyperfine structures have been observed in the 3s4s 3S, — 3s4p ®P0,1,2 transi- 
tions of PIV witha 21 ft. concave grating and in the 4s ?S1, — 4p*P%. transition of 
PIII with a transmission echelon. The half-widths of the lines, excited in an 
electrodeless discharge, were reduced by mixing neon with the phosphorus vapor 
and by using a condenser of small capacity. Their intensities were increased by 
discharging the condenser through the exciting coil nine times per half-cycle by 
means of a rotating spark gap. The structures of the P IV lines confirm J = 
1/2 and show that yp, the nuclear magnetic moment, is positive. w= 1.15 + 0.05 
n.m. and w = 1,11 + 0.1 n.m. are obtained from the P IV and P III structures, 
respectively. 












Introduction 





The nuclear magnetic moments of nearly all the nuclei lighter than Yt 
with odd atomic number and odd mass number are known. P*! is one excep- 
tion. Its spin, J = 3, has been determined from alternating intensities in 
the band spectrum of diatomic phosphorus vapor (1, 10). Its magnetic 
moment has been investigated by Tolansky (14): he found no resolvable 
hyperfine structure in any of the 13 lines of the second spectrum of phosphorus 
studied by him, and concluded that the nuclear gyromagnetic ratio, g, is 
anomalously small. However, these lines of P II would have very narrow 
structures even for a g-value corresponding to a magnetic moment of the 











order of one nuclear magneton (3). 






The fourth spectrum of phosphorus is more favorable for a determination 
of the magnetic moment. The transitions 354s 3S, — 3s4p *P§,,2 in P IV 
(2, 13) should have hyperfine structures of the order of seven times larger 
than the widest structure in the lines studied by Tolansky. Accordingly 
these transitions were investigated. An electrodeless discharge was used to 
excite the fourth spectrum, and the effects of source conditions on the half- 
widths of the lines were studied. The reduction effected in the half-widths 
made it possible to resolve hyperfine structures and determine the magnetic 










moment. 






Experimental 





The light source was an electrodeless discharge in a mixture of phosphorus 
vapor and neon in a quartz tube 75 cm. long and 2.5 cm. in diameter. The 
discharge was observed through a plane window at one end. An exciting 
coil 50 cm. long was wound tightly on mica wrapped around the central portion 
of the tube. Neon was circulated through the discharge tube. The con- 
struction and operation of the circulating and purifying system are described 
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in the book by Tolansky (15). Freedom from hydrogen and mercury in the 
discharge tube is essential since both appreciably quench the higher spectra 
of phosphorus. 

A schematic diagram of the exciting circuit is shown in Fig. 1. The oil 
condensers, C, combined in series as shown, had a capacity of 0.002 uf. The 
exciting coil, D, consisted of 118 turns of 1.5 mm. copper wire, and had an 


Q900000000 
000000000 


3. 1. Exciting circuit. 


inductance of 27.5 wh. It was found empirically that smaller values of the 
capacity resulted in sharper lines, but reduced the intensity appreciably. In 
order to use small capacity without loss of intensity, the method of multiple 
’ discharge per half-cycle was adopted. The condensers were discharged 
through the coil by the rotating spark-gap, S, which was synchronous with the 


applied voltage. The number of discharges per half-cvcle was determined 
by the number of fixed terminals per quadrant, and as many as nine were 
used. With the above value of the capacity the condensers charged up to 
approximately the same potential before each discharge. The charging 
current was supplied by a 1.5 kva., 220 : 75,000-v. transformer. Its input 
was controlled by a 1.5 kva. autotransformer and was adjusted to give a 
sparking voltage of approximately 20,000 v. 


The best compromise between intensity and half-width for P 1V lines was 
obtained with the electrical conditions described above and a pressure of 1 mm. 
of neon. Since the discharge tube was operated at room temperature, the 
partial pressure of the phosphorus vapor was that of yellow phosphorus at 
room temperature. The effects of source conditions were not investigated in 
sufficient detail to warrant much speculation on the mechanism of the excita- 
tion process. However, the empirical results indicate that the line widths of 
the higher spectra depend more upon the velocities acquired by ions in the 
exciting field than upon other factors, and thus show that data relevant to 
the theory of discharges in gases at low pressure can be obtained from the 
study of line contours. 

The spectrograph was a 21 ft. concave grating in a modified Eagle mounting 
with a plate holder 12 ft. long. It was in a constant temperature room 
separated from the source room except for a small opening through which the 
light entered. The 3350 A region was recorded in the fourth and fifth orders 
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simultaneously. In the fifth order, which is the more intense, the theoretical 


resolving power is 400,000. 

Ilford Zenith plates were used and developed in Ilford metolhydroquinone 
developer. A Corning red—purple Corex filter was placed before the slit to 
reduce the background and the overlapping of orders. Although the plates 
were grainy, traces obtained with a Leeds and Northrup microphotometer 
were relatively smooth because a sufficiently long portion of the line was 
scanned to average out fluctuations. The separations, as well as the inten- 
sities of the hyperfine structure components, were determined from the traces. 










An accurate calibration of the emulsion was obtained from the nonuniform 
intensity distribution along a spectral line due to the astigmatism of the 
concave grating, as first suggested by Dieke (5). Dieke’s method is applicable 
only if the slit is uniformly illuminated and the grating is uniformly reflecting 
along the length of a ruling. Since these conditions were not fulfilled, his 
method was modified to incorporate the principle of the ““Two-Line Method”’ 
described by Churchill (4). On a concave grating spectrogram the relative 
intensity distribution along the length of a spectral line is the same for all lines 
in a narrow spectral region. Thus the densities of two selected lines measured 
at equal distances from the centers of the lines represent a constant intensity 
ratio independent of the distance from the centers. By measuring corre- 
sponding densities at a number of distances from the centers of the lines the 
characteristic curve can be determined if the intensity ratio of the two lines is 
known. If the intensity ratio is not known, the characteristic curve can be 
plotted only on an arbitrary log J scale. 

In this investigation the intensity ratio of the transitions *S, — *P? (3364 A) 
and 3%, — *P9 (3371 A) was used to determine the log J scale uniquely. 
This ratio was first measured using a characteristic curve obtained from a 
stepped-slit exposure on another plate of the same emulsion number, and was 
within 10% of the ratio 3:1 given by multiplet theory for LS coupling. This 
is the order of agreement expected when the standard is not on the same plate 
as the spectrum. Thus it was deemed justifiable to use the theoretical ratio 
3:1 to determine the log J seale of the characteristic curve obtained by the 
Two Line Method. The intensity ratios for hyperfine structure components 
and for the lines of the fine-structure multiplet as determined with this 
calibration are consistent with theory to within 5% as shown in Fig. 2. 






















This consistency is appreciably better than that obtained with a stepped-slit 
calibration, particularly at low densities. This is not unexpected since the 
Two Line Method uses the lines themselves, and hence is independent of 
source fluctuations and variation of emulsion sensitivity over the plate. 
Further, the background corrections can be made with greater precision. 








Results 










The intensity contours of the three P 1V lines studied are given by the heavy 
lines in Fig. 2, (4), (6), (c). They were obtained from the microphotometer 
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traces of the best plate. Traces of seven other plates showed no appreci- 
able departure in form from those of the best plate. Relative intensities 
obtained from several of these plates agreed with those obtained from the 
best plate to within 10%, and half-widths agreed to within 2%. In our 
judgment the contours obtained from the best plate are reliable to better 
than 5%, and are probably more accurate than a simple average of the con- 
tours obtained from all the plates. 


2 





(a) (Co) 
3371.2 A 6 3347.8 A 
3S,-°Po 35,-°P? 








= 

< (b) 

= : 3364.5 A ‘ 
0 

Bn 38,-°P, 





The dashed contours in Fig. 2 represent the theoretical hyperfine structure 
components, assumed to have a Gaussian distribution. They are based on 
the theoretical patterns shown in Fig. 3, (a), (6), (c). These patterns are 
drawn for J = } and a positive nuclear magnetic moment. ‘The relative 
intensities are shown in Fig. 3 by the heights of the lines representing the 
components and were obtained from the formulae of Hill (9). The hyperfine 
structure separations given in terms of g in Fig. 3 are calculated as described 
below. The theoretical line contours shown in Fig. 2 as light continuous 
lines are the sums of the dashed contours of the components. Only two 
experimentally derived data were used in plotting the theoretical contours: 
0.44 cm.-! as the half-width of a component, deduced from 3348 A; and 
0.54 (5) cm.~! as the separation of the hyperfine structure states of *Si, 
obtained from the measured separation of the peaks (0.49 + 0.02 cm.~) in 
3371 A by correcting for the small overlapping of the components. The 
intensity unit for both the theoretical and the experimental contours is the 


intensity of the stronger peak of 3371 A. 
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The presence of just two hyperfine structure components in the pattern of 
3371 A with a measured intensity ratio of 1.98 proves uniquely that J = 3, 
the value obtained by Ashley (1) and Jenkins (10) from the molecular spectrum. 


ee 


3371.2 A 3364.5A 3347.8A 
(a) (b) (C) 


Fic. 3. Theoretical patterns. 


1 


The hyperfine structure splitting of the 4S, level for J = 4 is given by (7, 8) 


where 


3 
AGS) =5- mi 


5 (Gs. + Gu) «Xo em; (1) 


K is a relativistic correction and has the value 1.02 for an s 
electron of phosphorus, 


Qns = a’ ee “Mm. (2) 
n; 


R is the Rydberg constant, 
a is the Sommerfeld fine-structure constant, 


Z; is the inner effective nuclear charge and equals 15 for an s 
electron of phosphorus. 


Zo is the outer effective nuclear charge, 
n, is the effective principal quantum number, 


g is the ratio of the nuclear magnetic moment in Bohr magnetons 
to the nuclear angular momentum in Bohr units. 


Que = 232.8 Zog/n? cm . (3) 
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The values of Z2/n? are different for the two electrons of 3s4s °S,, and are 
not obtainable directly. Since the Rydberg formula for the term value gives 
only Z2/n?, Zo or nm. must be determined explicitly for each electron. For 
the 3s electron Z3/n? is evaluated by equating it to the term value, in Rydberg 
units, of the 3s4s configuration of P IV relative to the 4s configuration of P V. 
A value of n, to be used with Z3/n? is obtained as follows: n, of the 3s 
electron of P V evaluated from Robinson’s data (13), T = 524,463 cm.-! 
and Zo = 5, is 2.287. This is an approximate value of n, for the 3s electron 
in the 3s4s configuration of P IV. It is true that Zo = 5 is too large for the 
3s electron, by perhaps 10%, when the shielding by the 4s electron is taken 
into account; but the term value, T = 524,463 cm.~!, taken from the P V 
analysis is also too large by 12% compared to the energy required to remove 
the 3s electron from the 3s4s configuration of P 1V, namely T = 460,384 cm.~'. 
These two errors tend to compensate. Since Zo probably is between 5 and 
4.5, the value of ”,. obtained in this manner should not be in error by more 
than + 5%. Using nm, = 2.287 and T = 460,384 cm.—! for the 3s electron, 


Zi/n? = T/R.n, = 1.83; . (4) 
Substituting this value in (3), 
d3, = 427 gcm.“ . (5) 


For the 4s electron of the 3s4s configuration, Zp is taken as 4, then n, is 
3.061, as calculated from the term value, 187,423 cm.~', of 3S; relative to 
3s 2S of P V. Thus from (3) 


a, = 130g cm." . (6) 
Substitution of (5) and (6) in (1) gives 
A (@S,) = 425 gcm.— . (7) 


The experimentally determined value of A (3S,) is 0.54 (5) cm.~', and the 
hyperfine structure is normal. Thus by Equation (7), g = + 0.00128. 
However, Fermi and Segré (6) have shown that if the quantum defects for 
an s-electron term sequence are not constant the right hand side of Equation 
(3) must be multiplied by a correction factor. From Robinson's data for the 
ns terms of P III and P V and for the 3sus terms of P IV the correction 
factor for the 3s electron is 1.016 and for the 4s electron 1.025. These 
corrections increase the coefficient of g in Equation (7) by 2%, and give 
g = + 0.00125 + 0.00005. The nuclear magnetic moment, y, in nuclear 
magnetons then is + 1.15 + 0.05 n.m. 


The separations of the levels of the 3s4p configuration are consistent with 
Houston's equations and show that the coupling departs somewhat from LS. 
The wave function for the triplet level with J = 1, designated as {*P?}, can 
be expressed as a linear combination of the triplet and singlet wave functions 
in pure LS coupling, designated as [*P®] and [!P%] respectively. Using 
Robinson’s term values (13), 


{8P0} = 0.994 [3P9] — 0.106 ['P9]. (8) 
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Equation (8) shows that the departure from LS coupling is small, and that 
the 3:1 intensity ratio assumed for the two lines used to calibrate the emulsion 
is valid to within 1%. However, the departure is significant in the calculation 
of the hyperfine structure splitting of the *P9 level. 


For this calculation (3) it is necessary to express the wave function as a linear 
combination of the jj wave functions of the two levels with J = 1. Designating 
these as [1/2, 1/2] and [1/2, 3/2], to indicate their origin in the jj subgroups, 


{3P9} = 0.873 [1/2, 1/2] + 0.488 [1/2, 3/2]. (9) 


The coefficients in Equation (9) when substituted in the formulae of Breit 
and Wills (3) for intermediate coupling give the following relation for the 
splitting of the *P? level: 


A(@P%) = (0.327 dis + 1.09 asp), (10) 


where 
= 427 gcm.—! from Equation (5), 
Avg “ 
“Tc” 10.6 gcm.-. 
207 cm.—', the separation between the J = 0 and J = 2 
levels, 
; = 13, obtained from Av by the spin doublet formula, 
= as 


Relativistic corrections are included in the numerical coefficients of Equation 
(10). Thus, A@P9) = 228 gcm.-' = 228 X 0.00128 cm.-', = 0.29 cm.—!, 
using the uncorrected value of g obtained from the +S, splitting. 

An experimental value of A(*P) was obtained by comparing the observed 
contour of 3364 A with the contours predicted using A(3S,) = 0.54(5) cm.-! 
and a series of values of A(?P?).. The comparison gave a fairly unique value, 
ACP?) = 0.27 + 0.04 cm.-!. The contour plotted using this value, which 
is 0.02 cm.~! smaller than the calculated splitting, matched the observed 
contour very well, including the hump on the high-wave-number side. This is 
evident from Figs. 3, (0), 2, (b); a small decrease in the splitting of *P9 
increases the separation of the two central components and produces a hump 
on the side of the contour. However, the difference between the calculated 
and the experimental value of A(?P9) is within the assigned error which takes 
into consideration both the precision of the matching and the 4% error in the 
value of A('S;) used in the analysis of the contour. Although this error is 
relatively large because of the incomplete resolution, the contour is compatible 
only with J = 3. A similar analysis of the contour of 3348 A gave A(@P%) = 
0.39 + 0.11cm.—. Bycalculation (3), A@P8) = 282 g cm.-! = 0.36 cm.~!. 
Thus, the hyperfine structures of these two lines, as well as that of 3371 A, 
are consistent with J = 1/2 and w = +1.15 + 0.05 n.m. 

The third spectrum of phosphorus was also investigated with a transmission 
echelon. The transition 4s*S,, — 4p ?P% (4246 A) (11, 13) was resolved into 
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two components with a separation of 0.11 + 0.01 cm.-' However, only one 
good spectrogram was obtained before two plates of the echelon became 
“uncontacted.”” The resolved structure on this spectrogram gave A(4s?S,,) = 
0.12 + 0.01 cm.-' = K.as,,. Thus by Equation (3) with a Fermi-—Segré 
correction factor of 1.026, 4 = + 1.11 + 0.1 n.m. 

During the writing of this report a letter was published by Pound (12) in 
which he gave + 1.1314 + 0.0013 n.m. as the nuclear magnetic moment of 
P*! obtained by the magnetic resonance method. Our values agree with 
his within the error set by us on the basis of experimental accuracy alone. 
In addition our results show that the nuclear magnetic moment is positive 
and confirm J = 1/2. 
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THE GAMMA RAYS OF RADIUM (B + C)! 


By K. C. MANN AND M. J. OZEROFF 








Abstract 


An experiment has been performed to verify the gamma-ray energy spectrum 
of Radium (B + C) as reported by Latyshev and associates. The spectrum was 
obtained by measuring the energies of photoelectrons ejected from a’ thin lead 
radiator by means of a thin-lens magnetic beta-ray spectrometer. The spectrum 
obtained failed to reveal the presence of transitions other than those reported 
by C. D. Ellis. Comparative results are presented. 





















Introduction 







The gamma rays of RaB and RaC have been reported in detail by Ellis 
and coworkers (5) from measurements of the energies of internal conversion 
electrons and of photoelectrons ejected from high atomic number materials. 
Later, Latyshev and his associates (1, 2, 8, 9, 10) measured the energy distri- 
bution of positrons formed in the process of internal pair production as well 
as internal conversion electron energies. They were able to report the detec- 
tion and measurement of all the gamma rays found by Ellis above 2 moc? 
(1.02 Mev.) plus 10 new gamma rays, the comparative intensities of the latter 
being high enough for reasonably easy detection. It is strange that these 
gamma rays were not found by Ellis, in spite of the extreme care which charac- 
terized his work, and for this reason it seemed advisable to investigate again 
the gamma-ray spectrum of radium in an effort to eliminate these uncertainties. 














Apparatus and Method 





The source used in this investigation was a silver ‘‘needle”’ containing 10 mc. 
of radium in equilibrium with its disintegration products. Gamma rays from 
this source were studied in a thin-lens magnetic beta-ray spectrometer of the 
type described by Deutsch, Elliott, and Evans (3). For this purpose the 
gamma rays were allowed to eject into the spectrometer photoelectrons from 
a thin lead foil 50 mgm. per cm.’ in thickness. The momentum distribution 
of the photoelectrons was measured by varying the focusing current through 
the spectrometer’s lens coil. For detection of electrons, a bell-type Geiger- 
Miiller counter was used with a mica window 4 mgm. per cm* thick. 










The spectrometer was calibrated using the well known F-line of Thorium B. 
The momentum corresponding to this conversion line has been accurately 
determined (4, 11) and is 








Hp = 1385.6 gauss-cm. 






The thorium source was prepared by precipitation as a sulphide from a thorium 
nitrate solution on a thin backing of mica. When used in the spectrometer 









1 Manuscript received March 15, 1949. 
Contribution from the Department of Physics, University of British Columbia, Van- 
couver, B.C. 
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the conversion line was quite sharp and gave a calibration constant for the 
spectrometer of 486.0 gauss-cm. per ampere. 

Observations of photoelectron peaks indicate that the resolution of the 
A(Hp) 


i” is 4% at half-peak intensity. 


spectrometer, 


Results 


Fig. 1 shows the curve obtained with the spectrometer, using the radium 
source. The solid curve is the composite curve of photoelectron peaks and 
Compton electron background, the latter arising from a brass absorber used 
to eliminate all primary beta rays. The dotted curve represents the Compton 
electrons alone, taken with lead radiator foil removed. 
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Fic. 1. The gamma ray spectrum of Radium (B + C). 


The values of the energies in Mev. assigned to each photoelectron peak are 
the gamma-ray energies, calculated as the sum of the photoelectron energy 
plus the shell binding energy of lead. All peaks, with one exception, were 
assumed to be due to photoelectrons ejected from the K-shell, since this is 
most probable. For these, the K-shell binding energy used was 0.088 Mev. 
(taken from X-ray absorption limits (7) ). In the one exception, a smaller 
peak following the 0.607 Mev. gamma ray was found to be 72 kv. more 
energetic and so was identified as an L-shell peak of the same gamma ray. 
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Table I compares the energies and relative intensities of the gamma rays 
found in this investigation with those reported by Ellis and coworkers, and 
by Latyshev. The relative intensities have been calculated on a basis of 
peak height above background corrected for variation of photoelectric cross 
section, using published curves (6). 









TABLE I 













’ 


Lines marked ‘“‘i.c.’’ denote internal conversion measurements. 
Lines marked ‘‘p.e.”” denote photoelectric measurements. 
Lines marked ‘‘i.p.’’ denote internal pair measurements. 





























































Ellis and coworkers Latyshev and coworkers Present investigation 
Gamma-ray energy,| Gamma-ray energy, Relative Gamma-ray Relative 
Mev. Mev. intensity energy, Mev. intensity 
.0529 i.c. 
0.2406 i.c. 
0.2571 ic. 
0.2937 i.c. 0.291 . 
0.3499 i.c. 0.352 - 
0.4260 p.e. 

.4980 p.e. 0.496 14 
0.6067 i.c. 0.606 i.c. 0.607 74 
0.766 = ic. 0.766 i.c. 0.781 13 
0.933 i.e. 0.933 i.c. 

1.41: Ep. 
25420. ix. 1.120 i.c. 12 9 
1.21 ip. 23 24 
1.238 wc: Le 2oe ats 
1.379 ix. 1.370 i.c. 
1.39: ‘ip. 49 1.40 23 
1.414 ic. 1.414 icc. 
1.52 ip. 29 
1.62 i.p. 22 
1.69 i.p. 17 
ey ae 100 1:47 100 
1.761 icc. 1.761 i.c. 
1:32" 2p, 17 
2.09 ip. 15 
22193. ic. 2.20 i.p. + i.c. 41 2.21 2s 
2.40 ip. 21 2.40 * 




























* Not estimated. 





The 2.40 Mev. gamma ray reported in this paper does not appear as a 
photoelectron peak. It was assumed that a weak gamma ray must be 
responsible for the Compton recoil electron distribution of the upper energy 
end of the curve. Choosing the end point of this distribution by inspection 
to be 2.14 Mev., the gamma-ray energy responsible can easily be shown to 
be 2.40 Mev., with an uncertainty of the order of +3%. 








As used in this investigation, the spectrometer is somewhat insensitive to 
gamma rays of energy less than 0.250 Mev. because of the long electron path 
with its attendant seattering and the relatively thick counter window and lead 
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radiator foil. Hence the low energy rays reported by Ellis did not appear 
here. In addition, conversion coefficients strongly affect the comparison of 
results. For example, the 1.414 Mev. transition is known to be almost 
completely converted and, as a result, it does not appear at all in our measure- 
ments. 

One faint peak from a 0.450 Mev. gamma ray was found which has not been 
reported previously. Its intensity is so low that it is certainly near the limit 
of detection of the spectrometer. 

An inspection of the comparative results indicates a radium gamma-ray 
spectrum along the Ellis pattern. Little evidence of the extra gamma rays 
reported by Latyshev and coworkers was found. However, all photoelectron 
peaks above 1 Mev. are relatively weak compared to the admittedly high 
Compton background and it is possible that the latter may have masked the 
missing lines. 

Techniques are now being developed in this laboratory in an attempt to 
reduce the Compton background, which if successful should increase the 
sensitivity of the spectrometer to low intensity gamma rays. The authors 
hope to be able to make a further report of this research in the near future. 
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